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1.0 SUMMARY 


Flight testing of aircraft with natural laminar flow experiments in the last 10 years has brought this 
means of drag reduction closer to practical application on small to medium-sized transports. To 
assist in forming a database relating the length of natural laminar flow to various wing design 
parameters, NASA initiated the Variable Sweep Transition Flight Experiment (VSTFE) in 1983. 
This experiment involved flight testing different gloves on the F-14 aircraft and correlating the 
measured transition results with analytical disturbance growth calculations. The Boeing Company 
has been under contract to support NASA in several phases of the VSTFE. This report describes a 
new laminar boundary layer stability analysis computer system, called the Unified Stability System 
(USS), developed by Boeing, and documents USS analyses of 24 flight test cases from the VSTFE 
and earlier experiments. 

Intelligent design of laminar flow airplanes requires the ability to predict the effect of the selection 
of airplane design parameters such as sweep, airfoil, and flight condition on the extent of laminar 
flow. Since the current state of the art for transition does not permit rigorous predictions, a 
semiempirical approach is required. The most successful correlations use disturbance amplifica- 
tion factors given by stability theory as the correlating parameter. This approach, originated in the 
1950s, is quite straightforward for simple situations, but in modem applications there are many 
ways in which this method could be applied: following a particular disturbance along a streamline, 
or at some preferred direction, or some preferred wavelength or frequency, combinations of dif- 
ferent disturbances, or some other approach as yet untried. The determination of the best method 
requires a statistical approach wherein a battery of reliable data are analyzed in various ways to see 
which is the most successful. With previous stability methods, this has involved an impractically 
large amount of labor. The analysis of a single transition data point requires many stability theory 
runs for each of the several approaches under consideration. Previous stability codes have re- 
quired much human intervention in the form of initial guesses or in the selection of physically 
meaningful results from among a large number of solutions, some of which were spurious. The 
result is that previous transition correlations have not been exhaustively compared against all data. 
The purpose of the USS is to automate stability calculations, making statistical comparisons 
feasible . 

The USS was developed to analyze spatial laminar boundary layer stability over a wide range of 
disturbance frequencies and orientations. In doing so it would overcome the limitations of old 
stability methods, which analyzed only enough disturbance conditions to satisfy particular growth 
philosophies. In the USS, these growth philosophies are included only in the final disturbance 
growth integrations, not in the stability analyses. This allows new philosophies to be tried without 
recomputing stability characteristics. In addition to the automatic calculation of a wide range of 
stability characteristics, the USS is intended to be easier to use than previous stability analysis 
codes. 

The USS consists of eight separate computer codes that set up the run as desired by the user and 
carry out the boundary layer analysis, stability analyses, and integrations of disturbance growth as 
necessary. The laminar boundary layer is determined using a finite difference infinite sweep meth- 
od with taper corrections and can account for forward sweep. The boundary layer stability is 
investigated using techniques developed by Dr. L. M. Mack and involves numerical integration 
of the disturbance equations through the boundary layer with disturbance eigenvalues that must 
be iteratively determined so the boundary conditions can be met. After the disturbance character- 
istics are calculated over a large range of frequencies and orientations, the total growth of various 
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disturbances using different growth philosophies is calculated. Presently, Tollmien-Schlichting 
(TS) growth at constant wave angle, crossflow (CF) growth using “irrotationality,” and “maxi- 
mum amplification” growth can be calculated. 

Although the USS was intended to be easy to use, the wide range of boundary layer conditions 
that it may encounter made this difficult. The resulting system is reasonably expensive in comput- 
er time, and a new user should read the usage sections of this report carefully and completely 
before using the USS. The system creates many files that can be used to examine the details of 
laminar boundary layers and their associated stability with regard to various disturbances. 

Fifteen conditions from the F-14 VSTFE flight tests were analyzed by the USS to determine 
disturbance growth and correlate it to the experimental transition points. These crossflow and 
Tollmien-Schlichting disturbance growth N-factors at transition were combined with N-factors 
from earlier laminar flow flight experiments on the F-lll and 757, which were reanalyzed using 
the USS. This group of transition N-factors scatter over a wide range, which makes their use as a 
prediction method for future laminar flow wing design questionable. The sensitivity of transition 
N-factors to small wing pressure field changes can account for some of this scatter. Some research- 
ers consider the addition of curvature terms to the boundary layer stability equations to be impor- 
tant in establishing a transition criterion with less scatter. Adding curvature terms arising from 
surface and streamline curvature to the USS is straightforward, but curvature arising from the path 
of the disturbance followed on a surface is dependent on the growth philosophy chosen by the 
user. This would involve an iterative procedure. 

Despite the scatter in the transition N-factor correlation presented in this report, the wide range of 
disturbances examined by the USS makes it a useful tool for analyzing the nature of disturbance 
growth in laminar boundary layers. The USS is recommended for general use in understanding 
the stability of laminar boundary layers on wings at speeds up to about a sonic freestream. 
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2.0 INTRODUCTION 


Laminar flow investigations in the last 10 years have increased the optimism about obtaining 
significant areas of laminar flow on small- to medium-sized transports with little or no boundary 
layer suction. Intelligent design of laminar flow airplanes requires the ability to predict the effect 
of the selection of airplane design parameters such as sweep, airfoil, and flight condition on the 
extent of laminar flow. To significantly add to existing flight data regarding laminar flow on 
sweptwings, NASA initiated the Variable Sweep Transition Flight Experiment (VSTFE) . The 
VSTFE is an effort involving the NASA Langley Research Center and the NASA Ames/Dryden 
Flight Research facility, with The Boeing Company providing design and analysis support under 
contract to Langley. Starting in 1986, flight tests were conducted at Dryden with an F-14 fitted 
with full-span upper surface wing gloves. Transition was measured for a wide range of flight condi- 
tions and wing sweep angles. Two gloves of different designs were examined: a “cleanup” or 
minimum modification glove and a 0.7-Mach glove designed by NASA Langley. Early results 
from flight tests with the cleanup glove are reported in reference 1. Boundary layer stability calcu- 
lations were done on the VSTFE cleanup glove flight data by Boeing and were combined with 
measured transition locations to provide a transition prediction method. Reference 2 reports on 
the Boeing support of the VSTFE carried out before the flight testing. Reference 3 describes 
Langley’s design of the 0.7-Mach glove. 

The application of these data to airplane design requires adjustment for differences between con- 
ditions of the experiments and those of the intended application. Ideally, this adjustment would 
be made through a rigorous theory. Since a rigorous theory does not exist, an empirical or semi- 
empirical approach is required. The most successful approaches so far developed use disturbance 
amplification factors given by stability theory as the correlating parameter. This approach, origi- 
nated in the 1950s, is quite straightforward for simple situations, but in modem applications there 
are many ways in which this method could be applied: following a particular disturbance along a 
streamline, or at some preferred direction, or some preferred wavelength or frequency, combina- 
tions of different disturbances, or some other approach as yet untried. The determination of the 
best method requires a statistical approach wherein a battery of reliable data are analyzed in 
various ways to see which is the most successful. With previous stability methods, this has involved 
an unpractically large amount of labor. The analysis of a single transition data point requires many 
stability theory runs for each of the several approaches under consideration. Previous stability 
codes have required much human intervention in the form of initial guesses or in the selection of 
physically meaningful results from among a large number of solutions, some of which were spuri- 
ous. The result is that previous transition correlations have not been exhaustively compared 
against all available data. 


2.1 OBJECTIVE AND SCOPE 

As part of the flight data analysis task of the VSTFE, NASA contracted with Boeing to improve 
and expand the capability of using linear stability theory to determine disturbance growth in swept- 
wing laminar boundary layers. The desirability of expanding on available methods was due in part 
to the different philosophies used by different researchers in the industry. 

One philosophy used presently (refs. 2 and 4) investigates two classes of disturbances, those more 
or less aligned with the local external flow (called Tollmien-Schlichting (TS) disturbances), and 
those nearly perpendicular to the local external flow (called crossflow (CF) disturbances) . This 
philosophy involves choosing a wave angle at which to analyze the TS disturbances. This is similar 
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to the “Option 6” growth calculation in the COSAL stability analysis program, reference 5. The 
angle is usually chosen to give the greatest growth throughout the range of important frequencies. 
In addition, this philosophy considers the zero frequency (stationary) crossflow disturbances to be 
the most important in causing transition and calculates the growth of stationary crossflow waves for 
which the component of spanwise wavelength is constant. This is the “irrotational” method de- 
scribed by Mack in reference 6. 

Another philosophy of calculating disturbance growth does not distinguish between TS and CF 
wave classes, but calculates the growth of disturbances of different frequencies at whatever wave 
angle gives the maximum growth rate at each point on the wing. This is the philosophy used in 
“Option 1” of COSAL, reference 5. As shown in figure 1, both philosophies just described require 
only a partial knowledge of the boundary layer stability characteristics. 

The objective of the stability analysis procedure described in this report is to automatically deter- 
mine stability characteristics over a wide enough range of disturbance orientation and frequency 
so either of the disturbance growth philosophies mentioned above or any new approach could be 
evaluated. 

The stability analysis system developed during this task was to be used in examining the new 
high-quality flight transition data taken during the VSTFE. Disturbance growth would then be 
calculated by the philosophies described above to determine their viability as transition criteria. 

2.2 APPROACH 

The laminar boundary layer stability analysis procedure described in this report consists of eight 
computer codes and is called the Unified Stability System (USS). Figure 2 shows the sequence of 

Accumulated disturbance growth determined by- 
a Constant wave angle, various frequencies (for 
Tollmien-Schlichting) 

b Zero frequency, various wave numbers, 

“irrotational” (crossflow) 
c Maximum growth rate, various frequencies 
and wave angles 

Typical boundary layer Wave 

stability characteristics ^ 

surface at one point 
on a wing 

b 




= 90 deg ana unique 


wavelength. 

Figure 1. — Philosophies of Determining Disturbance Growth 


Disturbance 
growth rate 



Note: For each ijr 

and a) there is 
an associated 
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Figure 2. — The Unified Stability System 

execution of the eight codes. A master program sets up the job control statements to carry out the 
calculations desired by the user. Three programs set up the input for the boundary layer analysis, 
carry out the analysis, and prepare the boundary layer information for stability analysis. The 
boundary layer analysis is a finite difference method that can account for conventional or inverse 
wing taper. The boundary layer analysis and stability analysis codes were based on existing codes 
and only modified as necessary to satisfy the requirements of the USS. 

Two different computer codes are used to calculate the boundary layer stability. The mathemati- 
cal basis and solution procedure in both are almost identical, but one is tailored to analyze distur- 
bances with wave angles below 70 deg and the other disturbances with wave angles between 72 deg 
and about 91 deg. This range of wave angles covers the disturbances of primary interest to stability 
analysts for most cases on sweptwings. In choosing a stability analysis code for the USS, two 
existing methods were examined: COSAL, a temporal technique using matrix solution methods 
(ref. 5) and MACK, a method that can solve for temporal or spatial stability using numerical 
integration through the boundary layer (ref. 7). After considerable work with both methods, the 
MACK method was chosen for calculating the spatial stability required in the USS. 

The mathematical development of the linear stability theory used in the MACK method parallels 
that for the Orr-Sommerfeld equation, but compressibility is included and the spanwise dimension 
is added. Disturbances in the boundary layer are characterized as having some wave length, \, 
frequency, to, wave direction relative to the local external flow, \|r, and a spatial growth rate. 
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dN/ds. N is the exponent of “e” in describing disturbance growth as: disturbance amplitude at 
some point = e N times (amplitude of that disturbance when it first starts to be amplified) . The 
introduction of that disturbance into the compressible, 3D boundary layer equations results in an 
eighth-order system of equations with four unknowns; X, o>, \|f, dN/ds. Typically, wave length or 
frequency and wave orientation are chosen by the user, and the other two unknowns are solved 
for. This is an eigenvalue problem because solutions exist only for certain combinations of the 
unknowns. In the MACK method, guesses for the eigenvalues start an iterative solution process 
that uses the Newton-Raphson procedure to refine the eigenvalues until the system of equations is 
solved to within some adequate tolerance. In the USS, the initial guessing of values that will 
successfully converge to eigenvalues has been automated based on the work of previous authors 
and experience with the MACK method. 

When spatial disturbance growth rates for different disturbance frequencies and orientations have 
been calculated along the wing surface, the integration to find total disturbance growth is done by 
two different programs in the USS. Both constant wave angle (TS) and “irrotational” (CF) distur- 
bance growths are calculated, and the “maximum amplification” philosophy is also implemented. 
The constant wave angle integration is similar to “Option 6" of COSAL, and the maximum 
amplification” integration is similar to “Option 1” of COSAL. However, with the USS, all the 
stability analyses are done before the disturbance growth integrations, so calculating disturbance 
growth for a range of frequencies does not require different stability program runs, as it does with 
COSAL. 

The programs of the USS generate numerous files. Some only transfer information between pro- 
grams, but several are available to the user for detailed examination and plotting of the laminar 
boundary layer and its stability characteristics. 

Data for 15 flight conditions from the VSTFE cleanup glove testing were used to analyze laminar 
boundary layer disturbance characteristics for different wing sweep angles, freestream conditions, 
and wing pressure fields. The USS was used to calculate disturbance growth rates over a wide 
range of disturbance orientations and frequencies. Then downstream amplification of the distur- 
bances using the constant wave angle approach for TS waves, the irrotational stationary wave 
approach for crossflow waves, and for some cases the maximum amplification approach was car- 
ried out. The amplification of TS and CF disturbances at the flight-measured transition point was 
plotted on a graph of Njs versus Ncf for each case analyzed in order to establish transition 
criteria that could be used for laminar flow wing design. To further define the transition criteria 
using a consistent analysis method, several cases from earlier laminar flow flight tests on an F-l 1 1 
(ref. 8) and a 757 (ref. 9) were reanalyzed with the USS. 
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3.0 SYMBOLS AND ABBREVIATIONS 


Abbreviations: 

CF 

TS 

USS 

VSTFE 

Symbols: 

c 

Cp 

H 



RecF 


RerR 




Re^* 

s 

u 

X 

z 

a 






crossflow. 

Tollmien-Schlichting. 

Unified Stability System. 

Variable Sweep Transition Flight Experiment. 


wing section chord length, 
wing section pressure coefficient. 

boundary layer shape parameter, 8*/0. However, H (crossflow) is 
defined for the crossflow velocity profile as the ratio of (1) height at 
which the crossflow velocity is greatest to (2) the height at which the 
crossflow velocity is Xo of its maximum value. 

freestream Mach number 

accumulated growth of a disturbance in the laminar boundary layer. 
Also called N-factor. 

Reynolds number based on maximum crossflow velocity, edge 
kinematic viscosity, and distance from the surface to the height at 
which the crossflow velocity is Xo the magnitude of the maximum 
crossflow velocity. 

Reynolds number based on freestream conditions and the length of 
laminar flow at a wing section. 

Reynolds number based on local boundary layer thickness, 8, edge 
velocity, and edge kinematic viscosity. 

similar to Re<> with 8 replaced by displacement thickness, 8*. 
distance along an airfoil section surface from the attachment line, 
a velocity component in the w x” coordinate direction, 
distance along an airfoil section chord line, 
distance perpendicular to the airfoil chord line, 
angle of attack. 

real part of the nondimensional disturbance wave number in the u x” 
direction, see figure 9. In the rotated coordinate system, this is also 
the magnitude of the real part of the wave number. 

imaginary part of the nondimensional disturbance wave number in 
the M x” direction, see figure 9. In a coordinate system with x aligned 
with the group velocity, (*i is the negative of spatial disturbance 
growth rate. 
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<*r 5 


fir 


£(F-S) 

8 

8 * 

0 

.A. 

v 

(P v )wall 

(PU)oo 

(pQ)e 

* 

* 


0) 

X 

Subscripts: 

CF 

e 

LE 

TR 

TS 


the spanwise component of a r . This is used only with the unrotated 
coordinate system, and therefore, equals /J r in that system. 

real part of the nondimensional disturbance wave number in the “z” 
direction. Equals zero by definition in a coordinate system with x 
aligned with the group velocity. 

Falkner-Skan boundary layer pressure gradient parameter. 

boundary layer thickness, defined to be the distance from the 
surface at which the local velocity is 99% of the local edge velocity, 

boundary layer displacement thickness. 

boundary layer momentum thickness. 

sweep angle. 

local kinematic viscosity. 

suction mass flow rate. 

mass flow rate of the freestream. 

local mass flow rate at the edge of the boundary layer, 
an angle that defines the disturbance orientation, see figure 9. 

an angle that defines the orientation of the group velocity of a 
disturbance. This is the direction in which the disturbance grows in 
space. 

disturbance frequency, 
disturbance wavelength. 


denotes a parameter associated with crossflow disturbances 

denotes the local value of a parameter at the edge of the boundary 
layer. 

denotes the value of a parameter at the wing leading edge, 
denotes value of a parameter at boundary layer transition. 

denotes a parameter associated with Tollmien-Schlichting (low wave 
angle) disturbances. 


Superscript: 


denotes a dimensional quantity. 
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The following symbols are only applicable to section 4.2. They are presented separately because 
section 4.2 has the most complicated set of symbols in this report. 


Symbols: 


AjO) 


A 0 


C v» Cp 


D 


H (crossflow) 


L 

M 

P 

P 

R 

r cf 


r <5* 


S 

T 

t 

U (U.V.W) 


components of the characteristic vector of the compressible stability 
equations outside the boundary layer. The subscript denotes the 
component of the solution vector and the superscript denotes the 
number of the independent solution associated with this component. 

Initial disturbance amplitude in equation (21). 

Components of the eight first-order differential equations that the 
compressible stability equations can be reduced to. The subscripts 
vary from 1 to 8 and denote the position of a in the 8 by 8 matrix 
of terms. 

specific heat at constant volume and pressure, respectively. 

— , D 2 is the second derivative. 

dy 

components of the rate of strain tensor, see equation (2a). 

crossflow boundary layer shape factor. The ratio of (1) the height at 
which the crossflow velocity is greatest to (2) the height at which the 
crossflow velocity is Xo °f its maximum value. 

the imaginary number \ . 

the complex wave number vector. Without the subscript, it refers to 
the real part of the wave number. 

reference length. 

Mach number, 
pressure, mean value, 
pressure. 

Reynolds number. 

Reynolds number based on <5 10 , maximum crossflow velocity, and 
edge kinematic viscosity. 

Reynolds number based on 6 *, edge velocity, and edge kinematic 
viscosity. 

density, time fluctuating part. 

a “normalized” redefinition of Z done by the Gram-Schmidt ortho- 
normalization algorithm, see equations (18) and (19). 

temperature. 

time. 

the three components of mean velocity. U can also be used as the 
mean velocity vector. 
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u (u.v.w) 


Xj (x,y,z) 


Zi 


a 

P 

P(F-S) 

7 

A( ) 

5* 

810 

Sij 

0 

0 

K 

X 


P 

a 


the three fluctuating components of velocity, u can also be used as 
the fluctuating velocity vector and, when used with the bar (-), is 
the total time-varying velocity vector. 

the three component directions, x can also be any of the three 
orthogonal directions when used with a subscript. 

redefined variables in the stability equations that form a system of 8 
first-order differential equations. The index varies from 1 to 8. 

component of the complex disturbance wave number in the x 
direction. 

component of the complex disturbance wave number in the z 
direction. 

Falkner-Skan boundary layer pressure gradient parameter, 
ratio of specific heats. 

denotes a finite change in the parameter in the parentheses, 
boundary layer displacement thickness. 

height in the boundary layer at which the crossflow velocity is Xo the 
magnitude of the maximum crossflow velocity. 

Dirac delta function. 

in the normal mode form of disturbances, © is the function in the 
exponent of e, see equations (7) and (21). 

temperature, time-varying part. 

coefficient of thermal conductivity. 

second coefficient of viscosity. 

characteristic values that form the solutions to the system of 
first-order differential equations of stability outside the boundary 
layer. The index varies from 1 to 8 although only i = 1,3,5 and 7 
are physically reasonable. 

coefficient of viscosity. 

density. 

Prandtl number. 


Tjj stress tensor. 

an angle that defines the disturbance orientation, see figure 9. 

\p an angle that defines the orientation of the group velocity of a 

disturbance. This is the direction in which the disturbance grows in 
space. 

a) disturbance frequency. 
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Subscripts: 


CRITICAL 


i.j 


i 

r 

1 


value of a parameter at which the disturbance it is associated with 
has its maximum growth rate. This also applies to the value of a 
parameter when there is some disturbance that is unstable. 

both i and j are used to denote one of the three component 
directions. 

the imaginary part of a parameter, 
the real part of a parameter. 

denotes value of parameter at the edge of the boundary layer. 


Superscripts: 

* denotes a dimensional quantity except for 8*. 

denotes the fluctuating part of a variable. 


Special designators: 


caret, denotes the part of a variable that is a function of y only in 
the normal mode disturbance designation. See equation (7). 

overbar, denotes a parameter that varies with time. 

tilde, denotes a change of variable definition in equation (8). 
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4.0 USS THEORETICAL DEVELOPMENT 


The three main functions of the Unified Stability System (USS) are laminar boundary layer analy- 
sis, boundary layer stability analysis, and disturbance growth calculations. The first two of these 
functions use computer programs that were developed previously and described in the literature. 
The boundary layer method is described in detail in reference 10, and therefore only a brief 
outline of it is repeated here. The stability analysis uses a method developed by Dr. L. M. Mack, 
and it is explained in reference 7. However, that reference gives a more general background of 
stability theory than is used in the USS, so the development of the method that pertains to the USS 
is presented in this report. The disturbance growth calculations use computer code developed 
completely under this contract and are explained herein. 

4.1 BOUNDARY LAYER ANALYSIS 

The boundary layer program is essentially the same code, with some improvements and correc- 
tions, that was delivered to NASA in 1979 and reported in reference 10. It is a general-purpose 
3D finite-difference code that has been applied to wings, fuselages, inlets, wind tunnel walls, and 
propfan blades. In the present application, the code is used in the infinite-span-wing mode with 
taper effects, in which it solves the special form of the 3D boundary layer equations applicable to 
the coordinate system shown in figure 3. It is assumed that the isobars of the sweptwing flow lie 
along lines of constant percent chord, or that pressure is constant along the spanwise x-coordinate 
lines. 

It is further assumed that the boundary layer solution is similar along x-coordinate lines and that 
the boundary layer thickness varies with the local chord in a prescribed way (vt~ for laminar 
flow). The x-direction derivatives in the 3D equations can then be expressed in terms of deriva- 
tives normal to the surface, and spanwise finite differences are not required. The 3D solution can 
be generated by integration along a single chordwise arc. 

The same similarity arguments are applied to the flow in the neighborhood of the leading edge 
attachment line, leading to a special form of the attachment line boundary layer equations that the 
program solves to start the solution near the leading edge. Thus, for laminar flow the equations 
solved are equivalent to those solved by the special-purpose code of Kaups and Cebeci (ref. 1 1) . 
Running the 3D code in this mode requires specialized inputs. The z-coordinate is defined as arc 
length measured from the attachment line along the intersection of the wing surface and a sphere 
of radius r 0 centered at point A, normalized by c, the local chord at the intersection of the arc and 
the trailing edge. The metric coefficients consistent with this choice are as follows: 

h, = 1 
h 3 = c 
K,a = 0 
K 31 = — l/r 0 

The boundary conditions required along the arc are velocity vectors at the boundary layer edge 
that are consistent with the tapered wing similarity assumption. These are calculated by the input 
preprocessing program, BLGL, using airfoil section geometry and a pressure (Cp) distribution 
from a 2D or 3D outer flow solution or from experimental measurements. Because the attachment 
line of the flow generally falls between input data points, a special curve-fitting technique is used in 
the preprocessor to locate the attachment line and define flow quantities in its neighborhood, so 
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the boundary layer code computes correct flow derivatives and obtains a correct solution for the 
attachment line flow. For most laminar flow applications it is not adequate simply to assume that 
the attachment line coincides with the leading edge. 

Subsequent stability analysis of the boundary layer requires that the velocity profiles be defined 
along a different chordwise path than the arc along which the boundary layer code generates the 
solution for the velocity field. A postprocessing program, LAMSD, performs the transformation 
of the velocity profiles (a simple normal coordinate stretching) to the spanwise locations desired 
for the stability analysis. 
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4.2 STABILITY ANALYSIS 


The material presented in this section is an abbreviated form of what Dr. L. M. Mack gave in his 
AGARD paper (ref. 7). The important mathematical points pertinent to the USS stability analysis 
codes MKMOD3 and MKMOD5 are outlined. For more detail, reference 7 should be consulted. 

Compressible stability theory starts with the derivation of the governing equations from the Navier- 
Stokes equations for a viscous, heat-conducting, perfect gas, which in dimensional form are: 
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— r+ — r 
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r ij= 2 r*ij + |(X--r)f k *k-^ v 


(2a) 

(2b) 


The asterisks denote dimensional quantities, overbars denote time-dependent quantities, and ten- 
sor summation convention has been adopted. The equations are, respectively, of momentum, 
continuity, energy, and state. 

The stability equations are obtained from the Navier-Stokes equations by the following procedure. 
First, all quantities are divided into mean flow and fluctuating terms. With primes used to denote 
fluctuations of the transport coefficients, 


£T* = U* +U*, J P'=/ > * +/>*, 

T* = T* + 0*, p~* = p* + r*, (3) 

£*=**+*'*, X*=A*+A'*. 

where the first variable on each RHS is a steady mean-flow quantity, and the second is an un- 
steady fluctuation. 

Next, the equations are linearized, the mean-flow terms are subtracted out, and, finally, the paral- 
lel-flow assumption is made. The resulting equations are then made dimensionless with respect to 
the local freestream velocity U* 1f a reference length L*, and the freestream values of all state 
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variables (including the pressure). Both viscosity coefficients are referenced to M | . and k* is 
referenced to c p /i j, where c p is the specific heat at constant pressure. The transport coefficients 
are functions only of temperature so their fluctuations can be written 


, . dp \ , die 

p - — 0, k - — 

^ 1 dT i \dT 


) 0, X' = (— ) 9. 

I \ dT ) 


(4) 


Therefore, p., k, and X in the following equations along with p are mean-flow quantities, not 
fluctuations. The dimensionless, linearized x-momentum equation is 
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The y-momentum equation is 
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The z-momentum equation is 
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The continuity equation is 
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The energy equation is 
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The equation of state is 


p = r/Q + 6/T. 


(50 


d Si 

New terms that appear in these equations are M, y, and a = K , the Prantl number, which is a 
function of temperature. 

Equations (5) are the compressible counterparts of the incompressible stability equations and are 
valid for a 3D disturbance in a 3D mean flow. It should be noted that unlike most compressible 
stability analyses, equation (5e), the energy equation, is valid for a variable Prandtl number. The 
constant Prandtl number form is recovered by replacing k with ^ in the three terms in which it 
occurs. 

The boundary conditions at y = 0 are 

u( 0)=0, v (0) = 0, w(0)=0, 6(0) =0. (6a) 


The boundary conditions on the velocity fluctuations are the usual no-slip conditions, and the 
boundary condition on the temperature fluctuation is suitable for a gas flowing over a solid wall. 
For almost any frequency, it is not possible for the wall to do other than to remain at its mean 
temperature. The only exception is for a stationary, or near-stationary, crossflow disturbance 
when 6(0) = 0 is replaced by D0(O) = 0. The boundary conditions at y -* « are 

w(y), v 0). w(y), p(y), 6(y) are bounded as y -► » (6b) 

This boundary condition is less restrictive than requiring all disturbances to be zero at infinity, but, 
in supersonic flow, waves may propagate to infinity, and we wish to include those that do so with 
constant amplitude. 

Normal-Mode Equations. Except for the asymptotic suction boundary layer, most boundary 
layers grow in the downstream direction, and even for a wave of constant frequency, the distur- 
bances are all functions of x (and z in a general 3D boundary layer). What we have to deal with is 
a problem of wave propagation in a nonuniform medium. Because the complete linearized equa- 
tions are not separable, they do not have the normal modes as solutions. The most straightforward 
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approach is to simply set the nonparallel terms to zero on the grounds that the boundary layer 
growth is small over a wavelength, and it is the local boundary layer profile that will determine the 
local wave motion. This approach, called the quasi-parallel or locally parallel theory, has been 
almost universally adopted. It retains the parallel-flow normal modes as local solutions but is, of 
course, an extra approximation beyond linearization and leaves open the question of how impor- 
tant the admittedly slow growth of the boundary layer really is* We now specialize the disturbances 
to normal modes, 

[u,\,w,p,r,d\ T = [u(y), v(y), w(y),p(y),f{y), 0(y)] r exp[/ ^ J adx + fidz - wt j ], (7) 

where we have adopted the quasi-parallel form of the complex phase function. The normal modes 
may grow either temporally or spatially or both depending on whether o), or It, or both are com- 
plex. If a and (3 are real and w is complex, the amplitude will change with time; if a and (3 are 
complex and u> is real, the amplitude will change with x. The former case is referred to as the 
temporal amplification theory; the latter as the spatial amplification theory. If all three quantities 
are complex, the disturbance will grow in space and time. The original, and for many years the 
only, form of the theory was the temporal theory. However, in a steady mean flow the amplitude 
of a normal mode is independent of time and changes only with distance. The spatial theory that 
was introduced by Gaster (refs. 12, 13, and 14), gives this amplitude change in a more direct 
manner than does the temporal theory. 


When equations (7) are substituted into equations (5) and linear combinations of the x and z 
momentum equations formed for the variables, 

au=au+/Jw, dvv = aw-/3ti, (8) 


we obtain a system of equations. The momentum equation in the direction parallel to the wave 
number vector It is 


p [i(aU + pW - (o)au + (aDU + £DW)v ] = - i(a 2 +£ 2 ) {fi/yM x 2 ) (9a) 

A 1 

+ — [aD 2 u + (a 2 +0 2 ) (iDv- 2au)] + — — ~ (a 2 +fi i ){iD v-au) 

R 3 a 

+ -J ^{aD 2 U + PD 2 W)d+l ^-Dd + ^DTO \(aDU + pDW) 

R dT \ dT dT I 

+ —DT[aDu + /(a 2 + /? 2 ) v]|. 
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The y momentum equation is 


ip(aU + pW-a))v = - DP/yM 2 + ^-[2D 1 v + iaDii - (a 2 +/? 2 )v] 
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(9b) 


+ (D 2 v + idu) ( aDU + pDW)6 + 2 ^ DTD v 
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The momentum equation in the direction normal to k is 


p [i(aU + fiW - a>)aw + {aDW-pDU)\) = £ [d£> 2 vv - (a 2 +/3 2 )dvv] 
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The continuity equation is 


i(aU + PW - w)r + q{Dv + iau) + Opv = 0. 


(9d) 


The energy equation is 


p [i(aU + pW-w)d+ DTv ] = - (y- l)(Dv + /dti) 
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oR I tc 
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a 2 +0 


-[(aDU + pDW)aDu + {aDW-pDU)aDw) 


The equation of state is 


p = f /p+ 0/7. 


(9e) 


(9f) 


To reiterate, in these equations the eigenfunctions of the fluctuations are functions only of y and 
are denoted by a caret or a tilde; the mean-flow velocities U and W are also functions of y as are 
the other mean-flow quantities: density p(= 1/T), temperature T, viscosity coefficients p and X., 
thermal conductivity coefficient k, and Prandtl number. The specific heats are constant. The 
reference velocity for U and W is the same as for R and Mi, and the reference length for y is the 
same as in R. 
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Numerical Solution Procedure 

Eighth-Order System. Equations (9) are the basic equations of the compressible stability theory, 
but are not yet in a form suitable for numerical computation. For this purpose we need a system of 
first-order differential equations with the dependent variables defined by 

Z\ — clu + /?vv, Zj = DZ\ t Z 3 =v , (10) 

Z 4 =yS/yA/?, Z 5 = 0, Z 6 = DO, 

Zi = aw-f}u, Z 8 = DZ 7 , 

Equations (9) can be written as eight first-order differential equations 

8 

DZiiy) = ^ a ij(y) z i(y), (/= 1 . 8 ). ( 11 ) 

> 1 

and the fact that this reduction is possible proves that equations (9) constitute an eighth-order 
system. The lengthy equations for the matrix elements are listed in reference 7. 

The boundary conditions at the wall are 

Zi(0)=0, Z 3 (0) = 0, Z 5 (0)=0, Z 7 (0) = 0, (12) 

Z,(y), Z 3 (y ), Z 5 (y), Z 7 (y) bounded as y -> 

Uniform Mean-Flow Boundary Conditions. To solve equations (11) we also need to know the 
“boundary” conditions outside the boundary layer. In the flow outside the boundary layer, U = 
U lf W = W 1# T = 1, = 1, k = 1 fo u all y derivatives of mean-flow quantities are zero, and 

equations (11) reduce to a system of equations with constant coefficients. In spite of the greater 
complexity of these equations compared to those for incompressible flow, we are still able to arrive 
at analytical solutions. The lengthy derivation is given in reference 7. The exact freestream solu- 
tions are the ones to use to calculate the outer boundary values for a numerical integration of 
equations (11), but they do not lend themselves to a ready physical interpretation. For this pur- 
pose we examine the limit of large Reynolds numbers. The characteristic values simplify to 


A, .a ==F [a 2 + p - M\{aU x +p\V x -w) 2 ) x/2 , 

(13a) 

^3,4=^ [iRiaUi+pWx-a))]'' 2 , 

(13b) 

As , 6 =T [ioR{aUi+pWi ~u)V /2 , 

(13c) 

A?, 8 = A 3 , 4- 

(13d) 


We can now identify our solutions as, in order, the inviscid solution, the first viscous velocity 
solution, a viscous temperature solution that does not appear in the incompressible theory, and 
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the second viscous velocity solution. We shall only use the upper signs in what follows, as these are 


the solutions that enter the eigenvalue problem. 

The components of the characteristic vector of the inviscid solution are 

A 1 (1 > = -i(a 2 + /® 2 ) 1/a . ( 14a ) 

A 3 (1) = [a 2 + /3 2 -M?(al/, + /W, - a>) 2 ] 1/2 /(a 2 +/? 2 ) 1/2 (14b) 

A 4 (1) =i(aUi +pW x -co)/(a 2 + /? 2 ) 1/2 , (14c) 

Aj (1) =i(y- + /9W, - <u)/(a 2 + /? 2 ) 1/2 . (14d) 


The normalization has been changed to correspond to incompressible solutions. It can be noted 
that these expressions are correct when we set M, = 0. 

The components of the characteristic vector corresponding to the viscous velocity solution are 


> 

II 

(15a) 

A 3 (3) •-i/[iR{aUi +)3W, -w)] 1/2 . 

(15b) 

a 4 ' 3) = o ( a 5 (3) = 0. 

(15c) 


This solution is identical to the \ 3 incompressible solution only in the limit of large Reynolds 
numbers. 

The components of the characteristic vector corresponding to the viscous temperature solution 
are 


A,< 5) = 0, (16a) 

A 3 < 5 > -w) 1/2 /(iaR) l/1 . (16b) 

A 4 < 5 > = 0 . Aj (5) = 1. (16c) 

The components of the characteristic vector corresponding to the second viscous velocity solution 
are 

Ai ( 7) = 0 , A 3 (7) = 0 , A 4 <7) = 0 , A 5 (7) = 0 , (16a) 

A 7 < 7) = 1 . (16b) 

A 8 (7) = - [a 2 + £ 2 + iR(aUi + 0W\ ~ ( 1 6c) 


This solution is exact and is the same spanwise viscous wave solution as in incompressible flow. 
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We may observe that the viscous velocity solutions have only fluctuations of velocity, not of pres- 
sure or temperature. The velocity fluctuations in the x, z plane are in the direction of Ic for the first 
solution and are normal tolc for the second solution, which is periodic only in time. The viscous 
temperature solution has no velocity fluctuations in the x, z plane, or pressure fluctuations. We 
may regard these solutions as the responses to sources of u, w, and 6 ; and to emphasize this fact 
the respective solutions have been normalized to make these quantities unity. The second viscous 
velocity solution still has the interpretation of a normal vorticity wave, as in incompressible flow, 
but this wave cannot exist as a pure mode in the boundary layer (Squire mode) because of the 
dissipation term that couples the latter two of equations (11) to the first six equations. 

Since the early 1960s, the asymptotic theories developed by Tollmien (ref. 15) and Lin (ref. 16) 
have been largely superseded as a means of predicting numerical results in favor of direct solutions 
of the governing differential equations on a digital computer. The numerical methods that have 
been employed fall roughly into three categories: (1) finite-difference methods, used first by Tho- 
mas (ref. 17) in his pioneering numerical work on plane Poiseuille flow, and later, by Kurtz (ref. 
18), Osborne (ref. 19), and Jordinson (ref. 20), among others; (2) spectral methods, used first by 
Gallagher and Mercier (ref. 21) for Couette flow with Chandrasekhar and Reid functions, and 
later improved by Orszag (ref. 22) with the use of Chebyshev polynomials; and (3) shooting 
methods, used first by Brown and Sayre (ref. 23). All of these methods have advantages and 
disadvantages that show up in specialized situations, but they are all equally able to do the routine 
eigenvalue computations required in transition-prediction calculations. However, a shooting 
method used by Mack is what will be described here. 

The Shooting Method. Various integrators have been used to implement the shooting method. 
Perhaps the most common is some form of the Runge-Kutta method, but the Adams-Moulton and 
Keller box method have also been used. One choice that has to be made is whether to use a fixed 
or variable step-size integrator. The latter is better in principle, but it adds to the computational 
overhead, and thus to the expense, and it may be as difficult to construct a proper error test and 
then choose the error limits as it is to select the proper fixed step size. The integration method 
used in the USS is a modified fourth-order Runge-Kutta scheme described in reference 24, which 
eliminates some of the function evaluations required by the classical fourth-order method while 
preserving accuracy. This is a change to what Dr. Mack used in his computer code and is helpful 
in reducing computer time, since a high percentage of the effort of finding eigenvalues involves 
repeated applications of this integration scheme. 

Gram-Schmidt Orthonormalization. The early applications of shooting methods suffered from 
the problem of parasitic error growth. This growth arises because of the presence of a rapidly 
growing solution in the direction of integration that is associated with the large characteristic value 
X 3 in the freestream, which the numerical round-off error will follow. The parasitic error eventual- 
ly completely contaminates the less rapidly growing solution associated with the characteristic 
value X ! in the freestream. The essential advance in coping with this problem, which had previous- 
ly limited numerical solutions to moderate Reynolds numbers, was made by Kaplan (ref. 25). The 
Kaplan method “purifies” the contaminated solution by filtering out the parasitic error whenever 
it becomes large enough to destroy the linear independence of the solutions. 

A widely used method, which was originally developed for systems of linear differential equations 
by Godunov (ref. 26) and Bellman and Kalaba (ref. 27) and applied to the boundary layer stabil- 
ity problem by Radbill and Van Driest (ref. 28) and Wazzan, Okamura, and Smith (ref. 29), is 
that of Gram-Schmidt orthonormalization. This method has the advantage that it is easier to 
generalize to higher order systems than is the Kaplan filtering technique. However, the geometri- 
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cal argument often adduced in its support, that this procedure preserves linear independence by 
keeping the solution vectors orthogonal, cannot be correct because the solution vector space does 
not have a metric. In such vector spaces, vectors are either parallel or non parallel. The concept of 
orthogonality does not exist. Instead, the orthonormalization method works on exactly the same 
basis as Kaplan filtering: the “small” solution is replaced by a linear combination of the “small” 
and “large” solutions, which is itself constrained to be “small.” 

For the simplest case of a 2D wave in a 2D boundary layer, there are two solutions, Z (1) and Z (z \ 
each consisting of four components. In the freestream, Z (1) is the inviscid and Z (3) the viscous 
solution. Although this identification is lost in the boundary layer, Z (3) continues to grow more 
rapidly with decreasing y than does Z (1) . The parasitic error will follow Z (3) , and when the differ- 
ence in the “magnitudes” of Z C3) and Z (1) as defined by an arbitrarily assigned metric becomes 
sufficiently large, Z (l) will no longer be independent of Z (3) . Well before this occurs, the Gram- 
Schmidt orthonormalization algorithm is applied. The “large” solution Z (3) is normalized compo- 
nent by component to give the new solution 

5 (3) =Z (3) /{z (3r Z (3) } 1/2 , (18) 

where an asterisk refers to a complex conjugate and {} to a scalar product. The metric adopted 
for the vector space is the usual Euclidian norm. The scalar product of Z (1) and S (3) is used to form 
the vector 

5 (1) = [z (1)s -(5 (3), Z <1) )5 <3) ]/(^ 1) ^ 1) ) 1/2 , (19) 

which replaces Z (1) and where ? refers to the quantity in the numerator. 

The numerical integration continues with S (1) and S (3) in place of Z (1) and Z (3) , and when in turn 
S (3) exceeds the set criterion of, say, 10 s with single precision arithmetic and a 36-bit computer 
word, the orthonormalization is repeated. With homogeneous boundary conditions at the wall, it 
makes no difference in the determination of the eigenvalues whether the Z (l) or S (o 

are used. A 

linear combination of the two solutions satisfies the u(0) = 0 boundary condition, but the v(0) = 0 
boundary condition will in general not be satisfied unless a,/?, and to are eigenvalues of the prob- 
lem. 

Although the orthonormalization procedure has no effect on the method of determining eigenva- 
lues, it does complicate the calculation of the eigenfunctions. The solution vectors of the numeri- 
cal integration are linear combinations of the original solution vectors Z (l) and Z (3 \ and it is 
necessary to “unravel” these combinations. 

Newton-Raphson Search Procedure. The Newton-Raphson method has been found to be satis- 
factory for obtaining the eigenvalues. The boundary condition on Z (1) [orS (1) ] is satisfied at the 
conclusion of each integration by a linear combination of the two solutions at y = 0. The guess 
value of one of the eigenvalues is perturbed by a small amount and the integration repeated. Then 
the other eigenvalue is perturbed and the integration repeated again. 

The corrections to the initial guesses of the eigenvalues, for example a r and cti , are obtained from 
the residual v(0) and the numerical (linear) approximations to the partial derivatives by 

[dv r (0)/da r ] Aa r -[dv r (0)/<3a,] Aa, = -v r (0) , (20) 

[av/(0)/aa r ] Aa r -[av/(0)/aai] Aa, = -v,(0). 
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The corrected a r and a* are used to start a new iteration, and the process continues until A a v and 

Aai have been reduced below a preset criterion. The values in this criterion are a compromise 
between solution accuracy and computing time. These values were examined during the develop- 
ment of the USS and are given in the descriptions of output from programs MKMOD3 and 
MKMOD5 in section 7.0. 

Starting Guesses for the Eigenvalue Search. The search procedure for eigenvalues requires 
initial guesses for the eigenvalues. Depending on the nature of the boundary layer, the guesses 
may have to be quite accurate for the iterative search procedure to converge to actual eigenvalues 
rather than diverge. In developing the USS, a substantial part of the effort was automating 
successful eigenvalue guessing. 

The MKMOD3 program examines low wave angle stability. In this program frequency and wave 
angle are chosen, the imaginary part of £ is assumed to be zero, and the remaining eigenvalues to 
be found are the real and imaginary parts of a. Reference 29 has been extremely helpful in 
estimating the stability characteristics of low wave angle disturbances. These estimates are used 
to determine (1) when to begin searching for unstable disturbances during boundary layer 
growth, (2) what frequency range will likely have unstable disturbances, and (3) the guess for a r 
to use in starting the eigenvalue search procedure. 

Stability calculations begin in the MKMOD3 program when boundary layer Rg* is 500 less than 
the critical Rg* given in reference 29. The critical Rg* is a function of pressure gradient and is 
illustrated in figure 4. The Rg* and £(f-s) are characteristics of the boundary layer profile which is 
perpendicular to the local isobars, since the reference 29 results are for Falkner-Skan boundary 
layers. The value of 500 subtracted from the critical Rg* was determined by experience with 



Falkner-Skan Pressure Gradient 
Parameter, 

0 (F-S) 

Figure 4. — Critical Reynolds Number for MKMOD3 Program 
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MKMOD3. This helps ensure a successful start to the calculations because unstable eigenvalues 
are easier to converge on. 

The frequency range over which to search for eigenvalues is also found from information pres- 
ented in reference 29. Again, it is related to Rs* and 3 (f-S) by- 

log Which ~ log w LO w = / ( [log Ra* ~ log Ra *criticai\ » P<f-S)) ■ 

The function f is shown graphically in figure 5. This frequency range is distributed around a critical 
frequency found from an empirically determined equation, which is 

log wcrjt “ -1.46 log Ra* our + 0-37 


The a, initial guess is another value empirically derived from reference 29 and experience with 
MKMOD3 and is: 


a, (initial) - a rcwT “ ~ 

a rcR[T is a function of R^ crit as given in figure 6. A a r is a function of (log R<j. - log R^c^x) 

and P(F-S) and is graphed in figure 7. This guess is for the lowest frequency in the array of w and \Jr 
at the first point in the boundary layer that is analyzed because this is the starting point of the 
whole eigenvalue search. 
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Figure 6. — Nondimensional Wave Number at Critical Frequency and f?e 6 - CRjT 


The initial guess for a* is less critical for convergence of the search procedure and is 0.005. As 
eigenvalues are found in these early stages of the calculations, their values are used to get guesses 
of a r and c*i for calculations at different u and \|r and downstream stations in the boundary layer. 
This logic obviously works best when the frequency and wave angle increments are small and when 
the boundary layer changes are small between stations. 

The MKMOD5 program examines high wave angle (crossflow) stability. In this program, a r and 
wave angle are chosen, the imaginary part of g is assumed to be zero, and the remaining eigenva- 
lues to be found are a> and a* . As with the low wave angle stability, a reference was helpful in 
determining when the boundary layer was likely to have crossflow instability (ref. 30). The rela- 
tionship between H (crossflow) and a critical crossflow Reynolds number is shown in figure 8. 

Stability calculations are done in MKMOD5 if Rcf is greater than Rcf crit + 10 and continue 
until Rcf is less than ^cfcrit ~ 15 • These ranges were determined from experience with the 

program. Additionally, plots of spatial amplification rate for stationary (w = 0) crossflow distur- 
bances shown in reference 30 indicate that for low crossflow Reynolds numbers the most ampli- 
fied wave numbers are near 1.4 for wave numbers nondimensionalized by Sio- This observation is 
also used to start MKMOD5 stability calculations. 
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When stability calculations begin in MKMOD5, the initial eigenvalue guesses are based on experi- 
ence with the program as well as on published information. Experience has shown that when the 
calculations start at (R C f * r cf crit + 10) and a r = 1.4, that a good guess for co is 10" 5 and for c*j 

is 0.02 when the wave angle is slightly less (~ 1.0 degrees) than a critical angle that is determined 
from the boundary layer crossflow profile. This critical angle is defined as the direction for which 
the crossflow velocity profile has its inflection point at the point where the crossflow velocity is 
zero. In many cases these conditions of a r , i r t and w give nearly the most unstable disturbance in 
the crossflow 0J/~90 deg) region. When the initial eigenvalue search is successful, the MKMOD5 
program works its way to other wave angles and values of a r in order to outline the unstable 
crossflow region. Further details of this are given in section 7.0 under the paragraphs that describe 
the MKMOD5 output on the OPXQQYZ file. When the first boundary layer station has been 
analyzed, MKMOD5 follows a path similar to that of MKMOD3 in guessing at eigenvalues at 
following stations based on eigenvalues found at the previous station. 

Sixth-Order System. Equations (11) can be solved by numerical techniques, but the fact that 
there are 16 real equations (one each for real and imaginary terms) and four independent solu- 
tions means that the computer time required to calculate an eigenvalue can be substantial. It is 
therefore important to know if it is possible to make use of a system of lesser order. We note that 
for a 2D wave in a 2D boundary layer, the system already is of only sixth order, as there can be no 
velocity component, either mean or fluctuating, in the z direction. Is there an exact reduction 
available from eighth to sixth order? The answer, unfortunately, as mentioned by Dunn and Lin 
(ref. 31) and explicitly demonstrated by Reshotko (ref. 32), is no. 

We may observe from the coefficient matrix of equation (11) that the only term that couples the 
first six equations to the last two is the coefficient that comes from the last term of the energy 
equation (9e) and is one of four dissipation terms. It is the product of the gradient of the mean 
velocity normal to k and the gradient of the fluctuation velocity in the same direction. It was 
proposed by Mack (ref. 33) to simply set this term equal to zero and use the resultant sixth-order 
system for the calculation of eigenvalues. The numerical evidence is that except near the critical 
Reynolds number this approximation gives amplification rates within a few percent of those ob- 
tained from the full eighth-order system and is most accurate at higher Mach numbers. 

Wave Propagation in a Growing Boundary Layer. We have already discussed some aspects of 
this problem earlier in this section, and we have chosen to use the quasi-parallel rather than the 
non-parallel theory. In the quasi-parallel theory, the normal-mode solutions are of the form 

M(x,y,r,r) =A 0 w(y;x) exp [/©(*, y,z, r)] , (21) 


with similar expressions for the other flow variables. The slowly varying amplitude A(x) of the 
non-parallel solution has been set equal to the constant A 0 , and 


©(jc.z.O = { X a(x)dx+P(x 1 )z-w(x 1 )t . (22) 

We have left 0 and w as functions of the slow scale x 1 to make it clear that fl0/flx = a , just as for 
strictly parallel flow. The eigenvalues a, 0, and u> are locally related in what is called the dispersion 
relation and the eigenfunction u(y;x) is also a slowly varying function of x. Consequently, at each 
x a different eigenvalue problem has to be solved because of the change in the boundary layer 
thickness, or velocity profiles, or, as is usually the case, both. The problem we must resolve is how 
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to “connect" the possible eigenvalues at each x so that they represent a continuous wave train 
propagating through the growing boundary layer. 

In a steady boundary layer, which is the only kind that we shall consider, the dimensional frequen- 
cy of a normal mode is constant. For a 2D wave in a 2D boundary layer, 0 = 0, and the complex 
wavenumber a in the spatial theory is obtained as an eigenvalue for the local boundary layer 
profiles. The only problem here is the relatively minor one of calculating the wave amplitude as a 
function of x from the amplification rate, and we will discuss this in section 4.3. 

“Irrotational” Wave Propagation. When the wave is oblique, 0 # 0, and it is not obvious how 
to proceed, since a is a function of 0, as well as of x. How do we choose 0 at each x? The answer is 
provided by the same procedure as used in conservative wave theory. When we differentiate 
equation (22) with respect to x (not x,) and z, we obtain 

de/dx=a,dQ/dz=p , (23a) 

grad © = k c , (23b) 

where lc c is the complex wavenumber vector. Thus it follows directly that 


curl k c = 0, 


(23c) 


andk c is irrotational. This condition is a generalization to a nonconservative system of the well- 
known result for the real wavenumber vector in conservative kinematic wave theory. 

In the boundary layers on an untapered wing, the mean flow is independent of z. Consequently, if 
we restrict ourselves to spatial waves of constant p at the initial x, they can be represented by a 
single normal mode, because the eigenvalue a will also be independent of z. Therefore, according 
to equation (23c) the sought after downstream condition on p is 

P = const . (24) 

One caution is that if the reference length L* is itself a function of x, as it will be if L* = 8* for 
example, the argument has to be slightly modified and equation (24) refers to p* rather than 
to p. The untapered assumption is used in the disturbance growth calculations of the USS for the 
“irrotational” method, even though the boundary layer has been calculated using “infinite taper” 
assumptions. 

It still remains to specify the initial value of p. Naturally occurring instability waves in a boundary 
layer will be a superposition of normal modes with a spectrum over both to and p that will depend 
on the particular origin of the waves. It is probably only in a controlled experiment with a suitable 
wavemaker that a single normal mode can be excited. For example, the vibrating ribbon first used 
by Schubauer and Skramstad (ref. 34) in their celebrated experiment excites a spatial 2D normal 
mode with the frequency of the ribbon. It is also possible to conceive of wavemakers that excite 
single oblique normal modes in boundary layers that are independent of z. Such normal modes 
will have an initial which matches that of the wavemaker, and, because the wave can grow only 
in x, the initial fi\ must be zero. These normal modes are well suited for use in stability calculations 
for the estimation of the location of transition. In the calculations, /3 r is assigned as a parameter, ft 
is zero, and equation (24) controls the downstream values of ft ■ Not only do these normal modes 
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represent physical waves that can be produced by a suitable wavemaker, but they are also conven- 
ient to use in all calculations of normal modes, such as transition prediction where we are inter- 
ested in the largest possible growth of any normal mode. In earlier work on 2D planar boundary 
layers, the angle \|/ was chosen as the parameter to hold constant, rather than /? r , as the wave 
propagates downstream. Although a t is nearly constant in such boundary layers, it changes 
enough so that the assumption of constant \Jr is not equivalent to equation (24) . In the work on 3D 
boundary layers, equation (24) is applied to the spanwise wavenumbers, but the direction of the 
spatial amplification rate is nearly parallel to the local potential flow, or more accurately, in the 
direction of the real part of the group-velocity angle. 

This was identified in reference 35 by an argument based on the Nayfeh-Padhye amplitude propa- 
gation equation. For spatial waves the real part of the complex group velocity has x and z compo- 
nents — , — respectively. The direction of the real part of the group velocity, rp , with respect 
da T dfi r 

to the local potential flow is therefore 



In practice, this can be done using the finite differences 


AO) \ 

A a T Jfi r held constant 

if there is adequate knowledge of the stability characteristics at this point in the boundary layer. 
This is exactly what the USS is intended to do, so the finite difference calculations necessary to 
find group velocity angles are done after the arrays of w, a r , and (3 r are found for a point in the 
boundary layer. This procedure loses accuracy if the arrays of w, ot r , and p r are not defined by a 
dense number of points, but the arrays generally requested by the USS user are acceptable to 
define the group velocity angle, which is usually less than 10 deg. 

4.3 DISTURBANCE GROWTH CALCULATION 


da) 

Wr 


a T held constant 


da) 

da r 


The two programs of the USS that integrate disturbance amplification rate, dN/ds, along the wing 
surface to find the N-factors use similar mathematical techniques. They differ mainly in the inte- 
gration philosophy, which determines the disturbance frequency, w, and orientation, \|r, at which 
to pick the growth rate. The important mathematical properties of the programs will be described 
first. 


The N-factors are calculated by numerically integrating the 
dal rule. s 


N = 


dN 

ds 


ds 


following equation using the trapezoi- 




S 0 is the point at which a given disturbance begins to amplify. Although the path that is followed 
may be considered as lying along the surface of a streamwise section of a wing, because the ampli- 
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fication rates have been calculated along a streamwise section, the integration distance, As, is 
influenced by geometry, local velocity, and the disturbance group velocity. The following sketch 
illustrates this characteristic of the USS. 
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The direction used to determine As in the integration involves the average of group velocity angles 
and external streamline directions at points i and i + 1. This simplified method is a compromise 
between calculating stability and integrating along the disturbance growth direction, which is not 
known a priori, and the simplest approach of integrating strictly along the streamwise section 
surface. For many sweptwing analyses the streamlines do not depart significantly from the stream- 
wise direction, and the group velocity is generally within 10 deg of the local velocity direction. 

The integration process uses a straight line interpolation (extrapolation if necessary) to locate the 
point at which a disturbance first begins to grow. If a disturbance goes into a stable region after 
experiencing growth, the integration continues until the N-factor becomes negative. It then stops 
until there is again amplification of the disturbance. 

The amplification rate, dN/ds, used in the integration, is found from the tables of dN/ds created 
by the stability analysis programs. These tables represent discrete points on the stability surface 

dN 

/(**) 

ds 

calculated for each station on the section. The w and t | t are determined by the integration philoso- 
phy described later, and the resulting dN/ds is found from interpolation of the tables. Various 
interpolations are necessary and for most of these a second-order Lagrangian technique is used. 
In certain situations, however, the stability surface may not be adequately defined to allow a 
second-order interpolation so first order is used for interpolating some of the variables. Limited 
extrapolation may be useful and is allowed in certain situations, but because it may give unrealistic 
results, a warning message is printed with the integration results if extrapolation was involved in 
that result. 

The purpose of integrating disturbance growth rates along the wing surface is to find the total 
disturbance growth that can be correlated with boundary-layer transition. However, there are an 
infinite number of different disturbances that may be present in the boundary layer so some 
rationale must be reached as to which are to be used in the integration. As this document is 
written, no universally accepted philosophy is available regarding the disturbance integration so 
two commonly used philosophies are included in N-factor integration by the USS. 

The first philosophy investigates two classes of disturbances: those more or less aligned with the 
local external flow (these are the classical TS waves) and those nearly perpendicular to the local 
external flow (called crossflow disturbances). The growth of TS disturbances is calculated for 
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different frequencies, all of which have the same wave angle, \(r. For each frequency and wave 
angle there is an associated, unique wavelength at each point of the wing surface, see figure 1. The 
growth is calculated for several wave angles, and the maximum growth for the frequency/wave 
angle combinations is considered to be representative of the physical situation. Unfortunately, for 
some boundary layers the wave angle that gives maximum growth does not occur at a wave angle 
that is low enough to be considered in the TS regime. For those cases the investigator must set 
some rather arbitrary maximum limit on the wave angles of disturbances that will be considered 
TS. For the crossflow class of disturbances, growth rates for different frequencies are again inte- 
grated as they proceed downstream. For the crossflow, however, disturbances at a constant wave 
angle are not followed. Instead, an “irrotational" consideration described in section 4.2 is used, 
which results in following disturbances that have a constant spanwise component of wave number, 
a * . This rationale comes from a physical consideration of the disturbance's movement down- 
stream and is explained in more detail in reference 6. Of the infinite number of frequencies that 
are present, any one or several can be considered in the integration. Many investigators consider 
the zero frequency (stationary) disturbance to be the one critical for transition, but this point is 
still debated. 

The second disturbance growth philosophy can be called the maximum amplification approach. 
There is no distinction between TS and CF disturbances in this approach. In it, a disturbance of a 
given frequency is followed downstream using the maximum amplification rate for that frequency, 
considering the complete range of wave angles at each streamwise station on the wing. This inte- 
gration is repeated for several different frequencies. A weakness of this approach lies in whether it 
is physically reasonable to follow a disturbance whose wave angle may change drastically as it 
moves downstream. 

The two integration philosophies are pictured graphically in figure 1. N-factors using each of these 
philosophies are calculated by the USS codes INTTS2 and INTCF2. For the maximum amplifica- 
tion method, results from INTCF2 should be used for cases with moderate crossflow. INTCF2 has 
the ability to search through both crossflow and TS disturbances, if both are available, in its 
determination of maximum growth rate, whereas the INTTS2 code can only search the low-wave 
angle (TS) disturbances. 
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5.0 USS PROGRAMS AND FILES 


Numerous files are created by the Unified Stability System (USS) . The file names are determined 
by the user-specified case identification code, the wing surface analyzed, which program the job 
execution begins with, and a run version identifier. In naming the files, the following convention is 
used: 

X can be either *‘U” or “L” depending on whether the run involves an upper or 
lower wing surface analysis. 

Y the program step (1, 2, 3, 4, 5, 6 or 7) where job execution begins for this run. 

STEP EXECUTION BEGINS WITH PROGRAM 

1 BLGL. 

2 A411. 

3 LAMSD. 

4 MKMOD3 (TS disturbance analyses). 

5 MKMOD5 (CF disturbance analyses). 

6 INTTS2. 

7 INTCF2. 

Z A version number (0, 1 9) assigned to certain files. 

QQ Case identification code. 

Program RDECK 

Purpose: This code creates a job deck for a USS run. 

Input: RDECK is run interactively. The name of its input file is at the user’s discretion. This 

input file contains seven blocked record sets, which are inputs to the seven programs 
in the system. The first record of this file contains input for RDECK. 

Output: When RDECK is run, it prints a message to the screen informing the user that a “lo- 

cal” job file has been created, which can be submitted as a batch job. The name of the 
local file is JSXQQYZ, where the letters X, QQ, Y, and Z are as described above. 

Program BLGL 

Purpose: This program sets up the input for the boundary layer calculation by A411. This in- 

cludes curve fitting the input C p - x/c distribution, locating the attachment line, defin- 
ing the boundary layer grid, and calculating appropriate velocity components at the 
edge of the boundary layer. 

Input: File IBLGL (tape 1). This file is the second record of the user’s input file (see the 

input description section). 

Output: File I411UQQ (tape 2) or I411LQQ (tape 3). Input decks for the boundary layer code 

A41 1, upper and lower surfaces, respectively. 
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File CREFUQQ (tape 22) or CREFLQQ (tape 23). These files contain geometry in- 
formation relating the A4 1 1 solution stations to streamwise x/c positions on the upper 
and lower surfaces. 

File OPXQQYZ (tape 6). Printed output from BLGL and other programs. 

Program A411 

Purpose: This is the finite element boundary layer analysis code. For this system it is restricted 

to la min ar boundary layers but can analyze both standard and inverse wing taper. 

Input: File I411XQQ (tape 5). Inputs for A411 created by BLGL previously. 

Output: File BLXQQ (tape 4). The BL-file is a binary file containing boundary layer solution 

parameters. 

File VPXQQ (tape 10). This is a binary file containing detailed velocity and tempera- 
ture profile data through the boundary layer at each A411 solution station. 

File OPXQQYZ (tape 6). Printed output from A411 is also placed on the OP-file as 
noted previously. 

Program LAMSD 

Purpose: This interface program reads A411 boundary layer solution files and then generates 

input decks for the stability calculations. In addition, LAMSD can create plotting files 
containing boundary layer information. 

Input: File IBLGL (tape 1). This file is used as input for LAMSD as well as BLGL. 

File CREFXQQ (tape 2). This is the geometry reference file created by BLGL. 

File BLXQQ (tape 4). One of the binary files created by A411. 

File ILAMSD (tape 5) . This file contains the first and third blocked record sets of the 
user’s input file. The primary use of this file is to supply the interface code with the 
locations at which boundary layer information is to be passed to the stability codes. 

File VPXQQ (tape 10). The binary file generated by A411 containing boundary layer 
profile data. 

Output: File VPGPXQQ (tape 3) . If desired, this file can be created, and it contains boundary 

layer profiles in plotting format at user-specified locations. 

File SIXQQ (tape 7). This file contains boundary layer parameters and profiles in 
binary form for use by the stability codes. The boundary layer information is at user 
specified locations. 
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File BLGPXQQ (tape 8) . This file contains a condensed set of boundary layer param- 
eters at each chordwise station analyzed by A411. The information is in plotting for- 
mat. 

File OPXQQYZ (tape 6). Printed output from LAMSD is also placed on the OP-file. 
Program MKMOD3 

Purpose: This is the modified MACK code, which calculates eigenvalues defining the boundary 

layer stability at wave angles less than 70 deg. 

Input: File IMACKTS (tape 3). This is record set 4 of the user’s input file. 

File SIXQQ (tape 12). This is the binary file containing the boundary layer informa- 
tion required for the stability calculations. 

Output: File OTSXQQZ (tape 4) . One of the two files from MKMOD3 containing information 

regarding the calculation of each eigenvalue. The OTS-file is the longer of the two files 
and is useful for diagnosing program failures. 

File OPXQQYZ (tape 6). The more brief of the two printed files from MKMOD3 is 
put onto the OP-file. 

File TSIXQQZ (tape 8). The TSI-file contains stability information needed to calcu- 
late the TS N-factors. 

Program MKMOD5 

Purpose: This is also a modified MACK code that calculates the boundary layer stability in the 

72- to 90-deg wave angle range. 

Input: File IMACKCF (tape 3). This file contains record set 5 of the user’s input file. 

File SIXQQ (tape 12). As with the MKMOD3 code, this file contains the boundary 
layer parameters and profile information. 

Output: File OCFXQQZ (tape 4) . One of the two files from MKMOD5 containing information 

regarding the calculation of each eigenvalue. The OCF-file is the longer of the two 
files, and is useful for diagnosing program failures. 

File OPXQQYZ (tape 6) . The more brief of the two printed files from MKMOD3 is 
put onto the OP-file. 

File CFIXQQZ (tape 8). The CFI-file contains stability information needed to calcu- 
late the CF N-factors. 
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Program INTTS2 

Purpose: This code uses the stability information from MKMOD3, and calculates the N-factors 

at lower wave angles (Tollmien-Schlichting) using both constant wave angle and maxi- 
mum amplification methods. 

Input: File ITSI (tape 5). This file contains the sixth record set of the user’s input file. 

File TSIXQQZ (tape 8). This file contains the stability information calculated by the 
MKMOD3 code. 

Output: File OTSIXQQ (tape 6). Printed output from INTTS2 is put only on this file, not on 

the OP- file. 

File GTSIXQQ (tape 10). Output in plotting format for plotting N TS versus x/c. 
Program INTCF2 

Purpose: This code uses the stability information from MKMOD5 and if necessary from 

MKMOD3 to compute the crossflow N-factors (irrotational method) and N-factors 
using the maximum amplification method. 

Input: File ICFI (tape 5). This file contains the seventh record set of the user’s input file. 

File TSIXQQZ (tape 8). This file is used in this program as well as INTTS2, because 
the maximum amplification method may require stability information at low wave 
angles. 

File CFIXQQZ (tape 9) . This file contains the stability information calculated by the 
MKMOD5 code. 

Output: File OCFIXQQ (tape 6) . Printed output from the INTCF2 code is put on the OCFI- 

file, not on the OP- file. 

File GCFIXQQ (tape 10). Output in plotting format for plotting N CT versus x/c. 

At the end of a run, two additional files are saved on the user’s account: one is the dayfile from 
the run, named LSXQQYZ, and the second is the output file OSXQQYZ. The output file OS- 
contains a copy of the job cards used to control this run and the user’s input file. In addition, if an 
abnormal termination occurs, then the contents of the OP-file will also appear on the OS-file. 
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The following table summarizes the information in this section. 


Program 

RDECK 

BLGL 


A411 


Input files required 

User input file (record 1). 

IBLGL user input file (record 2). 
This file is not made permanent 
on the user’s account. 

I411XQQ 


Output file s generated 

JSXQQYZ (“local submit” file) 

I411XQQ 

CREFXQQ 

OPXQQYZ 

BLXQQ (binary) 

VPXQQ (binary) 

OPXQQYZ 


LAMSD 


MKMOD3 


MKMOD5 


INTTS2 


INTCF2 


IBLGL (same as used for BLGL). 
ILAMSD (user input file 
records 1 and 3). This file is 
not made permanent on the user’s 
account. 

CREFXQQ 

BLXQQ 

VPXQQ 

IMACKTS (user input file record 4). 
This file is not made permanent 
on the user’s account. 

SIXQQ 

IMACKCF (user input file record 5). 
This file is not made permanent 
on the user’s account. 

SIXQQ 

ITSI (user input file record 6). 

This file is not made permanent 
on the user’s account. 

TSIXQQZ 

ICFI (user input file record 7). 

This file is not made permanent 
on the user’s account. 

TSIXQQZ (if available). 

CFIXQQZ 


SIXQQ (binary) 
BLGPXQQ 
OPXQQYZ 
VPGPXQQ (optional) 


OTSXQQZ 

OPXQQYZ 

TSIXQQZ 


OCFXQQZ 

OPXQQYZ 

CFIXQQZ 


OTSIXQQ 

GTSIXQQ 


OCFIXQQ 

GCFIXQQ 


JSXQQYZ 
(batch job 
submittal file) 


LSXQQYZ 

OSXQQYZ 
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6.0 USS INPUT DESCRIPTION 


The input deck to the Unified Stability System (USS) is divided into seven separate records. If a 
certain case does not require the input in one or more records, the end-of-record marks are still 
required. 

The first record of input is primarily used to set up the job control deck for the run. 

Card 1 A descriptive card that tells the user what is input on the next card. 


Card 2 (format A7, 3X, A7, 3X, A6, 4X, A7). 

USERN: User number. 

PASSWRD: Password. 

CHARGE: Charge number. 

LIB: Account where the USS codes are stored. 

Example: 


■ 

Q 

□ 

ri 

□ 

□ 

H 

□ 

□ 

|Q 

m 

E 

|Q 

OQ 

E 

EQ 

EH 

EQ 

m 

S3 

m 

F2 

m 

si 

a 

a 

m 

a 

a 

a 

Efl 

m 

a 

S3 


S3 

EB 

B 

S3 

□ 

31 

EB 

31 

B 

ED 

S3 

B 

ED 

31 

31 

si 

a 

a 

SI 

S3 

1 

I 

u 

b 

3 

b 

B 

0 

1 

■ 

1 

Q 

Q 

B 

B 

S3 

B 

D 

a 

■ 

■ 

D 

B 

0 

B 

D 

B 

1 

■ 

1 

1 

B 

B 

D 

B 

B 

151 


1 

1 



1 



1 

1 


1 

■ 





1 
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Card 3 A descriptive card that tells the user what is input on the next card. 


Card 4 (format A40). 

USER: Mailing information for the job. 

Example; 


1 

2 

3 

4 

0 

6 

7 

e 

9 

10 


12 

13 

14 

15 

16 

17 

IB 

19 

20 

21 

22 

23 

24 

25 

26 

27 

58 

29 

30 

¥ 

1 

N 

X 

§ 


H 

0 

L 

D 


F 

0 

R 

□ 

R 

0 

z 

E 

N 

D 

A 

A 

L 




□ 
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Card 5 A descriptive card that tells the user what is input on the next card. 

Card 6 (format A2, 8X f 3F10.0). 


ID: 

PRI: 

VERSION: 

TIME: 


Two characters identifying the case. 

Priority to give to this run. Not used for NASA. 

A version number assigned to the output files generated by the stability 

codes. Valid entries are 0., 1., 2. f 

The CPU time limit for this run in seconds. 


Example: 


l 

2 

3 

4 

nn 

€ 

7 

a 

9 

10 

11 

[72 

13 

14 

15jl6 

17 

18 

19|20 

21 

22 

23 

24 

25 

Hi 

26 

29 

30 

31 

32 

33 

34 

36 

36 


1 

39 

40 

A 

A 



- 


□ 

— 

L_ 







I 



1 


1 




L 

□ 





1 

5 

0 

0 



□ 





8-U80 106-1 03 

Card 7 A descriptive card that tells the user what is input on the next card. 


Card 8 (format 6F10.0). 

FIRST: The number of the first code in the USS to be executed first for this 

run. 

1. BLGL. 

2. A411. 
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3. LAMSD. 

0. None of the first three codes need to be run. 

END: The number of the last (of the first three codes) code to be executed 

for this run. 

TS: Number to determine whether or not the MKMOD3 code will be run 

during this run. (1. for yes, 0. for no). 

CF: Number to determine whether or not the MKMOD5 code will be run 

during this run. (1. means yes, 0. means no). 

TSI: Number to determine whether or not the TS integrator code, INTTS2, 

will be run during this run. (1. for yes, 0. for no). 

CFI: Number to determine whether or not the crossflow integrator code, 

INTCF2, will be run during this run. (1. for yes, 0. for no). 

Examples. 


D 

B 

Q 

□ 

Q 

□ 

K2 

□ 

□ 

KB 

KB 

KB 

KB 

KB 

KB 

KB 

KB 

KB 

KB 

S3 

m 

m 

m 

03 

m 

m 


m 


EE 

S3 

EE 

S3 

EB 

EB 

EB 

EB 

23 

EB 

S3 

31 

03 

a 

B 


S3 

CB 

03 

S3 

si 

SI 

a 

a 

a 

a 


■ 
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■ 

■ 

■ 

■ 

1 
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■ 

■ 

n 

■ 

■ 

■ 

■ 

■ 

■ 

1 

1 

i 
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■ 

1 

i 

i 

i 

i 

i 

1 

■ 

D 

1 

1 

1 

1 

1 

■ 

1 

1 

1 

D 

■ 
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1 

1 

1 

1 

D 

i 

1 

1 

T 

'h 

e 

al 

DC 

>v 

e 

e: 

<a 

nr 

‘P 

le 

i 

S 

a 

fi 

in 

r 

un of all codes in the system. 













E 

IT 

T 

4 

T 
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The above example is a run of only the MKMOD3 and INTTS2 codes. The first three codes must 
have been run previously to create the necessary input files for MKMOD3. 


Card 9: A description card that tells the user what is input on the next card. 

Card 10: (format A7, 3X, A7, 3X, F10.0). 

DFILE: The name of this input file. 

GP: The name of the A488 GP-file if C P - x/c input will come from an A488 

run. If the C p ’s are input in this file, put the word “NOTHING” in this 
space. 

SURF: Number to determine which wing surface is to be analyzed in this run. 

(1. for lower, 0. for upper). 

Example: 
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The second record of the input file is used by the BLGL code to set up inputs for the boundary 
layer calculation. 

Card 1: A descriptive card that tells the user what is input on the next card. 

Card 2: (format A80). 

TITLE: The title associated with this run. 

Example: 
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Card 3: 


A descriptive card outlining the input on the next card. 


Card 4: (format 4F10.0). 


EMINF: 

RBINF: 

TINFK: 

ALPHA: 

Example 


Freestream mach number. 

Chord Reynolds number in the freestream direction. 
Freestream static temperature in degrees Kelvin. 

Angle of attack. This parameter is not presently used. 


1 

[7 

3 

0 

5 

6 

7 

8 

0 

10 

1 1 

12 

13 

1* 

15 

0 

17 

18 

19 

20 

£1 


23 

24 

25 

s 

27 

28 

29 

30 


0 


8 







1 

0 

5 

0 

0 

0 

0 

0 




2 

1 

8 

_ 

J 



_ 






Card 5: A descriptive card outlining the input on the next card. 
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Card 6: (format 5F10.0). 


SWP1: 

XC1: 

SWP2: 


XC2: 

XCHORD: 


Sweep at the forward part of the wing in degrees. Must be greater than 0. 
x/c at which the above sweep is applicable. 

Sweep at the rear part of the wing in degrees. Cannot be exactly equal 
to SWP1 but may be greater than SWP1, which indicates “inverse” taper 
(forward sweep). 

x/c at which the above sweep is applicable. 

Chord length in feet in the streamwise direction. 


Example.; 
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Card 7: A descriptive card. 


Card 8: (format 4F10.0). 


ETA2: 


ETAMAX: 


FETA: 

FSUC: 

Exampki 


A boundary layer grid parameter. This is the grid spacing at the wall in 
decimal fraction of displacement thickness. The value 0.012 is recom- 
mended. 

A boundary layer grid parameter. This value is the maximum distance 
from the wall to which the boundary layer will be calculated in units of 
displacement thickness. A value of 8.0 is recommended. 

The number of points defining the boundary layer profiles. A value of 
70. is recommended. Maximum value is 100. 

The number of points that will be used to define the suction distribu- 
tion. The value 0. implies no suction. Less than or equal to 300. 
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If there is suction, the next set of cards are input. 


Card 8 A: A descriptive card. 

Card 8B: (format 2F10.0). 
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SSUC: s/c points defining the suction distribution on the wing. These values 

should start at s/c = 0. and end with a value greater than the s/c expected 
at the trailing edge. This distribution is fit with straight lines between 
defining points. 

VSUC: Nondimensional suction rate. '.PAwaiiAP U A . A negative value de- 

notes suction. 

Example: 
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Card 9: A descriptive card. 

Card 10: (format F10.0). 

Put the value 2.1 on this card. 


Card 11: 

Card 12: 
FNPTI: 
SWPCPI: 

FINPTC: 

FNSEC: 

Example: 


A descriptive card. 

(format 4F10.0). 

The number of points that will be used to describe the C p - x/c distri- 
bution. Should be less than 250. 

The correction factor sweep in degrees that will be applied to the C p 
and z/c input. This is used to apply simple sweep theory (i.e., if a “nor- 
mal” pressure distribution is available) the BLGL program can do the 
adjusting of C p by cos 2 -A. and z/c by cos Jl . 

Number indicating where the C p information will come from. 

0. - C p - x/c input follows. 

1. - C p - x/c are input from an A488 GP-file. (Not available) 

Number indicating the A488 wing section to be used if the 

Cp - x/c input is to come from an A488 GP-file. (Not available) 
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Card 13: This is a descriptive card required if C p - x/c are input in this file. 

Card 14: (format A40). 

FORM: The format that will be used to read in the x/c, z/c, and C p values. The 

format is enclosed in parentheses. 
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Example: 
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Card 15: This is a descriptive card, required if the C p ’s are input in this file. 


Card 16: (format as specified on card 14). 


XCI: x/c, z/c, C p input as sets on the streamwise airfoil. This set must start 

ZCI: at the trailing edge lower surface and go completely around the section 

CPI: to the trailing edge upper surface. Maximum of 250 sets of values. 

Example: 
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The third record of the user input file contains information used by the interface program, 
LAMSD, to provide boundary layer information to the stability programs at the wing stations 
designated by the user. 

Card 1: A descriptive card outlining what is input on the following card. 


Card 2: (format 3F10.0). 


FLAG: 


0 . 

1 . 

2 . 

3. 


FSAV: 


FGGP: 


A number that determines how the boundary layer profiles from A4 1 1 
will be chosen for input to the stability codes. 

Choose the first N stations. N is specified in the following field. 
Choose all stations generated by A411. 

Choose N profiles from stations that are closest to the s/c values specified 
by the user later in this record. This is the recommended method. 
Choose N stations selected by a geometric spacing criterion that is denser 
near the leading edge. 

The number of stations for which boundary layer profiles will be saved. 
If FLAG = 0. and FSAV = 0., the first 50 stations generated by A41 1 
will be saved. However, only 40 stations can be calculated in one run of 
the stability codes. 

The number of stations where profiles will be saved on a plotting file. A 
value of 0. means no profiles are saved. The maximum number of sta- 
tions that can be saved are all those generated by A411. 
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Example: 
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Card 3: A descriptive card required if the first number on the previous card was 2. 

Card 4: (format A40). 

FORM: The format used to read in the s/c values where boundary layer profiles 

will be saved for the stability codes. The format is enclosed in parenthe- 
ses. This card is required only if the first number on card 2 is 2. 

Example; 
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Card 5: A description card required if the first number on card 2 was 2. 

Card 6: (format as specified on card 4). 

SCPUT: s/c values where boundary layer profiles and other parameters will be 

saved for use by the stability codes. 

Example.: 
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Card 7: 

Card 8: 

FORM: 

Example: 


8-U80 1 06- 1 1 6 

A descriptive card required if the number in field three of card 2 is not zero, 
(format A40). 

The format that will be used to read in the locations whose boundary 
layer profiles are to be put onto a plotting file. The format is enclosed 
in parentheses. 


( 


> 


I2|i3[l4h5h6|l7 
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Card 9: A descriptive card required if the number in field three of card 2 is not zero. 

Card 10: (format as prescribed on card 8). 

FGGPL: The numbers of the A411 spanlines (chordwise locations) at which the 

boundary layer profiles will be put into a plotting file. 
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Example; 
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The fourth record of the user input file contains the input for the MKMOD3 code. 


Card 1: (format 415, F10.0). 


LFRST: 


LLAST: 


NPSI: 

NFR: 

DELPSI: 


Example: 


A number that is the wing station to start checking the local Re,j*» and 
if it is beyond some threshold (which is a function of Falkner-Skan 0) 
the MKMOD3 code starts calculating stability information. This sta- 
tion number is referenced to those that were put on the Si-file, not all 
those produced by the A411 program. 

The number of the station to stop doing stability calculations. This sta- 
tion number is also referenced to those put in the Si-file. There can be 
no more than 40 stations analyzed by the stability codes in one run. 
The number of wave angles to calculate stability for at each station 4 < 
NPSI < 8. 

The number of frequencies, oj, to calculate stability for at each station. 
4 < NFR < 10. 

The wave angle increment in deg. Recommended values are between 
about 12 deg to 20 deg. This may have to be reduced for difficult cases 
such as adverse pressure gradients. 
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The rest of this record is made up of namelist input for MKMOD3. 
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Namelist 

PSIB3D: 


NUNIT: 

NPR3D: 

0 

1 

2 

NTABPR 

0 

1 

2 

CKOUT: 

0 


STD AT: 

Group velocity angle with respect to the local external flow. This effects 
the stability results only slightly and is usually near 0. The recommended 
input value is 0. and the default value is 0. 

Tape unit on which the longer portion of printed output and optional 
output will be placed. Use 4 only. 

Optional printing during eigenvalue iteration. 

Minimum print (default). 

First level of extra printing. 

Second level of extra printing. 

Optional printing of the boundary layer profiles. 

No profile printing. This is the default and recommended value. 

Print profiles in F-format. 

Print profiles in E-format. 

More optional printing. 

None. This is the default and recommended value. 
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1 Level 1 extra printing. 

2 Level 2 extra printing. 

Esamptei 

'X”^1T”1T7273 14 7? 16 17 18 10 20 21 22 23|24 2S 26 27 26 

i A T PS I B 3 D = 0 . , N U Nil T = 4 


Namelist SADAT: 
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NSOL: 


Example: 


The maximum number of iterations allowed in the eigenvalue searches. 
The default and recommended value is 9. 

This variable controls whether the dissipation terms are included in the 
solution. A value of 3 is recommended and results in the dissipation 
terms being ignored. If NSOL is not input, dissipation will be included. 


2|3|4|S[6|7|6|0 1iO|i1|12|13|14|iS1i6|17|Tb 19 20 21 22 23 24 2S 2 

$ ADAT N I = 9 , N S 0 L = 3 $ 


27 12812® 130 131 132113 k 
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Namelist STAB: 


The number of evenly spaced intervals into which the boundary layer 
profiles are divided for integration in the MKMOD3 code. The default 
value is 100. The minimum number of steps that does not significantly 
effect the accuracy of the results is about 80. The maximum allowed is 
200. However, the user is cautioned that a case that involves difficult 
eigenvalue convergence may be assisted somewhat by using more inte- 
gration steps (i.e., increasing Nl). This is the only variable in TAB that 
may require changing at present. 


Example: 


| 9 |l0|llh2ll3|l4|lSi16[l7jia|i9|20| 


TAB N 1*80 
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The fifth record of the user input file has the input for the MKMOD5 program. 


All the namelist input described for the MKMOD3 input is identical for MKMOD5 . However, the 
first card is different. 


Card 1: (format 415, 3F10.0). 

LFRST: A number that is the station to start checking the local crossflow Rey- 

nolds number. If it is greater than some threshold (which is a function 
of the crossflow H parameter), the MKMOD5 code starts calculating 
stability information. This station number is referenced to those put on 
the Si-file. 

LLAST: These three quantities are defined in the MKMOD3 input. 

NPSI: 

NFR: 

ALPH AI : This is the guessed value for the nondimensional disturbance growth rate 

eigenvalue to be used in starting the stability calculations. For a run start- 
ing from the leading edge, the recommended value is 0.02. If the case 
being run might start at a station that is only slightly unstable, a smaller 
value will be required. For many cases the accuracy of this guess is not 
critical. 


44 




FRINP: This is the guessed value for the nondimensional disturbance fre- 

quency eigenvalue to be used in starting the stability calculations. This 
value is related to the following input parameter, DPSIST. For cases 
that start at the leading edge, a value of 0.00001 is recommended when 
DPSIST is near 2. If DPSIST is smaller, FRINP should be proportion- 
ally smaller, too. For cases that start at downstream locations, a smaller 
value of FRINP is appropriate. 

DPSIST; This is the incremental wave angle in degrees that is subtracted from 
the critical wave angle to define the wave angle at which stability calcu- 
lations will start. It is desirable to have a small (< 2.) value for DPSIST, 
and its value influences what is to be used for FRINP, as mentioned 
above. 1. to 2. is usually a good range of values to use for DPSIST. 


Example! 


1 

nn 

0 

4 

T\ 

[51 

7 

T 

0 

10 

1 1 

12 

13 

14 

15 

16 

17 

18 

1# 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 





1 




2_ 

2 





_6 





_6 


_0 


0_ 

2 





_ 


0 


0 

_0 

_0 

0 

1 






2 

_L 

0 











— 


8-U80106-123 

The sixth record of the user input file has input for the IN I I S 2 program. 


Card 1: (format 15). 


ITYPE: 


1 

0 


Card 2: 


A number that determines the integration philosophy to be used to find 
the N-factors. 

Integrate disturbances of different frequencies at the wave angle for 
maximum growth rate. 

Integrate disturbances of different frequencies at a user input wave 
angle. 

(format 15, 9F10.0). 


If card 1 has a 0, then this card is required. Otherwise, it isn't. 

NIPSI NIPSI’s the number of different wave angles for which N-factors will 

PSII: be calculated. The PSII are each of the wave angles in degrees. A 

maximum of nine wave angles can be used. 

This record is set up to do two integrations of N-factors; one with the maximum growth rate 
method and one with the constant wave angle method. 


Example; 
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The seventh record of the user input file has input for the INTCF2 program. 
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Card 1; (format 15). 

METHOD: The first field contains a number that determines the integration 

philosophy to be used to find the N-factors. 

1 Integrate disturbances of different frequencies at the wave angle for 

maximum growth rate. 
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0 Integrate disturbances of constant dimensional spanwise component of 

wave number and of a user-specified frequency (the “irrotational” 
method) . 

Card 2: (format 15, 9F10.0, /, 5X, 9F10.0). 

If card 1 has a 0, then this card is not required. 

NIFR: The number of different frequencies for which N-factors will 

be calculated. 

FRI: Each of the different frequencies in Hertz. A maximum of 

18 frequencies can be used. 

This record is set up to do two integrations of N-factors; one for the maximum amplification 
method and one with the irrotational method. 


Example: 
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7.0 USS OUTPUT DESCRIPTION 


This section describes the output of the various Unified Stability System (USS) programs. The file 
naming convention is defined in the Programs and Files section. 

RDECK Output File JSXQQYZ 

This file is made up of CDC job control statements for this particular run of the USS. The combi- 
nation of programs to run and files to use and create are reflected in this job stream. 

BLGL Output Files I411UQQ or I411LQQ 

One of these two files is created when the BLGL program is run, depending on whether the upper 
or lower surface of a wing is being analyzed. 

The first line on the 1411-file is the run title, and the second has flag settings for A411 which are 
fixed for the USS applications. The third line has these parameters, in order: 

1. The negative of the freestream Mach number. 

2. Chord Reynolds number. 

3. Freestream static temperature in degrees Kelvin. 

4. This fourth number does not apply to laminar boundary layer calculations. 

5. Prandtl number. 

6. This sixth number does not apply to laminar boundary layer calculations. 

7. The last two parameters are internal A411 flags and are not changed in the USS application. 
The fourth line has the following parameters, in order: 

1. ETA2 from the user input file. 

2. ETAMAX from the user input file. 

3. The rest of these parameters are fixed for the USS application. 

Line five has a fixed output control flag, and line six has the number of wing sections and number 
of chordwise divisions to be analyzed. There is only one section for the USS application, but the 
chordwise stations to be analyzed varies around 220. 

Line seven always has a 0. for the USS application, and line eight starts a set of input that gives the 
ZI coordinates of the chordwise stations to be analyzed. These values always start at zero. ZI is 
defined in figure 3. 

After the ZI input are two lines that have values not applicable for the USS. 
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The file is completed by a set of lines that are in groups of three. The values on the first two lines of 
these three are fixed for the present application. On the third line the first, second, sixth, and 
seventh parameters are variable. The first is the local external velocity component parallel to the 
sweep nondimensionalized by freestream velocity. The definition of the velocity components is 
illustrated in figure 3. The second number is the local external velocity component perpendicular 
to the sweep. The sixth value on this line is the curvature of the chordline used for the A411 
analysis. It would be zero for a nontapered case. The seventh value is surface suction defined as 
(pv) waI1 /(pQ) e . Note that the definition of this parameter is different from the user input suction 
parameter, VSUC. 

This three-line group of input is repeated for each chordwise station to be analyzed by A411. 
BLGL Output Files CREFUQQ or CREFLQQ 

One of these two files is created when the BLGL program is run, depending on whether the upper 
or lower surface of a wing is being analyzed. The first line has (in order) 

1. The total number of stations on the chord at which A411 may compute the boundary layer. 

2. The sweep in degrees at the wing trailing edge. If the case has inverse taper, this will be the 
conventional definition of wing sweep minus 180 deg. 

3. The last two numbers are not important for this program application. 

Following the first line are columns of information for each A411 solution station. The first col- 
umn is the number of the station. The second is the ZI coordinate as defined in figure 3. The third 
column is x/c, and the fourth is z/c. The next two columns are geometric information. Column five 
is the negative of the sine of local sweep and column six is the direction cosine in the planform 
plane of the surface direction vector along the section used by A4 1 1 ; see figure 3 . Column seven is 
the ratio of local streamwise chord to the user-specified streamwise chord; see figure 3. Column 
eight is the local sweep in degrees. If the case has inverse taper, this will be the conventional wing 
sweep -180 deg. The last column is s/c of each solution station but on the user-specified stream- 
wise section. 

BLGL Output on the OPXQQYZ File 

There are several groups of output from BLGL on the OP-file. Several are of minimal use and will 
not be explained here. The case title starts the output. The next important group has column 
headings that start with I, ZI, XC, etc. This group contains information on the upper surface of 
the user-supplied section. The “I” column contains the numbers of the A411 analysis stations. ZI 
is defined in figure 3 and is nondimensionalized by the chord of the user-supplied streamwise 
section. XC is x/c. XB, YB, and ZB are A411 coordinates defined in figure 3, nondimensional- 
ized like ZI. XL AM and ZLAM are angles in the A4U coordinate system as shown in figure 3. 
Note how these values are affected by inverse sweep. CAOCB is the ratio of “local” streamwise 
chord to the chord at the user-defined section. QQQ is the magnitude of the local external velocity 
nondimensionalized by the freestream velocity and CPQQQ is local Cp. The CPQQQ values are 
the result of curve-fitting done by BLGL of Cp versus s/c, so this dense array of points in this 
output group can be very useful for confirming the accuracy of the curve fit by plotting CPQQQ 
versus XC and spotting the user-supplied input points on the plot. The total number of stations 
input to A4 1 1 is around 220, and the density of them is greatest at the leading edge and at sudden 


48 



changes in C p . This output group is followed by several unimportant ones. The next useful one is 
the lower surface station information. Both upper and lower surface station information are out- 
put, even though only one of the surfaces will be analyzed per run by A411. 

The lower surface station output is followed again by some groups of unimportant output. The 
next noteworthy output is a repeat of the 141 1-file output, explained above. The last output from 
BLGL placed on the OP-file is a repeat of the CREF-file that was explained previously. 

A411 Output Files BLXQQ and VPXQQ 

These two files are in binary format and are only used to transfer information to the next program 
in the USS. 

A411 Output to the OPXQQYZ File 

The A411 program prints details of the boundary layer solution at each station, starting at the 
attachment line. It is important to note that the “UNIT RE" referred to in the output is actually 
chord Reynolds number. Most of the numerous other parameters whose values are printed are 
defined in reference 10. The exceptions are defined below: 

For the attachment line output only, spanline 1, there are several special parameters printed. 


They are: 


DSTR3Z 

The partial derivative of DSTR3, with respect to the surface 
coordinate perpendicular to the attachment line, ZI. 

HSSS 

The ratio DSTRS/THSS. 

TH13Z 

The partial derivative of TH13, with respect to ZI. 

CSTAR 

((H3) (UEB) 2 ]/[(kinematic viscosity at the edge of the bound- 
ary layer) (WEBZ)]. 

WEZBX 

0(0(WEB)/d(ZI))/0XI, XI is along the attachment line. 

CPZZ 

d 2 (pressure coefficient)/dZI 2 . 

K13Z 

The rate of change of spanline curvature. For the USS applica- 
tion, this is 0. 

For all spanlines (actually chordwise stations) after the attachment line, the additional parameters 
printed that are not defined in reference 10 are defined here. 

BETASMAX 

The maximum crossflow angle in the boundary layer. The defi- 
nition for BETAS is not conventional. See reference 10. 

BETASMIN 

The minimum crossflow angle in the boundary layer. 

DSTTOT 

Displacement thickness based on the local total velocity in the 
boundary layer. 
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THTOT 


OMEGABY, 

OMBSQRX, 

OMBSQRZ 


Momentum thickness based on the local total velocity in the 
boundary layer. 

These parameters only apply to applications that use a 
rotating reference frame, and therefore are 0.0 for the USS 
application. 


When A411 reaches laminar separation, a note will be printed at each station saying: SOLUTION 
FORBIDDEN AT THIS STATION. KFORB(I) = the station number at which separation oc- 
curred. This concludes the pertinent A411 output to the OP-file. 


LAMSD Output File VPGPXQQ 

This is an optional file that can be created if the user desires. The file is formatted to plot velocity 
profiles calculated by A411. The profile information from any of the A411 chordwise stations can 
be written onto this file as chosen by the user-supplied input to LAMSD. The profile parameters 
written onto the VPGP-file are— 

ETA Distance above the wing surface nondimensionalized by 8 of 

the velocity profile in the vertical plane which contains the ex- 
ternal velocity vector. 

U Velocity profile in the vertical plane which contains the exter- 

nal velocity vector, nondimensionalized by the magnitude of 
external velocity. 

UP d(U)/d(ETA). 

UPP d 2 (U)/d(ETA) 2 . 

W Velocity profile in the vertical plane perpendicular to the ex- 

ternal velocity vector. Positive is outboard and the magnitude is 
nondimensionalized by the magnitude of the external velocity. 

WP d(W)/d(ETA). 

WPP d 2 (W)/d(ETA) 2 . 

T Static temperature nondimensionalized by boundary layer 

edge temperature. 

TP d(T)/d(ETA). 

TPP d 2 (T)/d(ETA) 2 . 

The profiles at each station make up one set of output on the VPGP-file, and the number just after 
the KI printed at the beginning of each set is the A411 chordwise station number. 
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LAMSD Output File SIXQQ 

This file is in binary format and is used to transfer the boundary layer parameters and profiles to 
the MKMOD3 and MKMOD5 programs. User-specified input to LAMSD determines which sta- 
tions from A411 have information passed on to the stability programs. 


LAMSD Output File BLGPXQQ 

This is another file formatted for plotting, and it contains several of the more important boundary 
layer parameters for each station analyzed by A411. The information extends from the first sta- 
tion behind the attachment line to the location where A4 1 1 experienced laminar separation or to 
the trailing edge. 


The parameters printed to this file are defined below: 


ZI 


The chordwise surface coordinate in the A41 1 coordinate sys 
tern, see figure 3. 


XOC 

QEB, 

BETA1E, 

DSTR1, 

DSTR3, 

TH11, 

TH33 

RXFLO 


HI33 

HCF 

REDSTRZ 

BETAF 

CROSSM 


x/c of the station. 

These parameters are all defined in reference 10 or have been 
defined earlier in this section. 


Local crossflow Reynolds number. The length used is the 
greatest height in the boundary layer at which the crossflow 
velocity is Xo of its maximum value. The velocity value used is 
the magnitude of the maximum crossflow velocity. 

This is the incompressible boundary layer shape factor for the 
profile in the A411 “chordwise” direction. 

Crossflow shape factor. The height in the boundary layer 
where the crossflow velocity is greatest divided by the greatest 
height where the crossflow velocity is Xo * ts maximum value. 

Boundary layer Reynolds number based on displacement 
thickness. The profile used is the A411 “chordwise” one. 

Falkner-Skan beta parameter, calculated using the velocity 
profile in the A411 “chordwise” direction. 

Maximum crossflow velocity in the vertical plane perpendicu- 
lar to the external velocity vector. Positive is outboard. This is 
nondimensionalized by the magnitude of the external velocity 
vector. 


51 


LAMSD Output Placed on the OPXQQYZ-File 

The first printing that LAMSD does to the OP-file tells the user where A4 1 1 stopped its boundary 
layer calculations because of laminar separation. The second set of output starts with the user-spe- 
cified title for this run and whether the upper or lower surface is being analyzed. Following this, 
there is a summary of parameters supplied to the stability analysis programs. They are defined 
below. 


CHORD 

SWEEP 1 (DEG) 

SWEEP 2 (DEG) 

RADIUS 

FREESTREAM MACH 
FREESTREAM CHORD RE 


XCHORD from input. 

SWP1 from input. 

SWP2 from input. 

Radius in the A411 coordinate system; see figure 3. 
EMINF from input 
RBINF from input. 


STATION The number of the station, as picked by LAMSD from A4 1 1 

results to be transferred to the stability codes, is given here and 
also the total number of those stations. 

A411 STATION The number of this A411 station picked by LAMSD. 


XOC 


x/c of this station. 


zoc 

SOC CHORDWISE 


CHORD CORRECTION 
DEL THETA (RAD) 


LOCAL SWEEP(DEG) 

POINTS IN PROFILE 

DISPL THICKNESS 
(DSTRZ) 

B.L. THICKNESS 
(DELTA) 


z/c of this station. 

The s/c of this station measured from the attachment line on 
the user-input streamwise section. 

This is CA/CB; see figure 3. 

This is an incremental sweep angle in radians and referenced 
to the leading edge sweep. A negative number is given when 
inverse taper is input. 

Sweep at this station. 

This is FETA from the user-supplied A411 input. 

Displacement thickness, in feet, of the velocity profile in the 
vertical plane perpendicular to local sweep. 

Boundary layer thickness, in feet, to where the velocity is 
99% of the edge velocity. It is calculated using the profile 
in the vertical plane which is perpendicular to local sweep. 
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DELTA/DSTRZ 


Ratio of the above two parameters. 


DELTA/CHORD 
SHAPE FACTOR 

FALKNER-SKAN BETA 

A411 TURBULENCE 
FLAG 

FLOW ANGLE (DEG) 

PHI CRIT (DEG) 

ETA OF CROSS- 
FLOW MAX 

ETA WHERE CROSS FLOW 
IS 0.1 OF MAX 

CROSSFLOW REY 

REY DSTRZ 
MACK REY 

EDGE VEL (UE) 

EDGE TEMP (DEG R) 


Boundary layer thickness divided by airfoil chord. 

Boundary layer shape factor, 8*/0. It is calculated using the 
profile in the vertical plane which is perpendicular to local 
sweep. 

This is calculated using the shape factor just above. 

0. for the USS application. 


The angle in degrees of the local edge flow measured from the 
local chordwise direction, positive outboard. 

This is the so-called critical angle, the angle at which the sec- 
ond derivative of crossflow velocity with respect to the “verti- 
cal” coordinate is zero at the same vertical distance as the 
crossflow velocity itself is zero. It is referenced to the plane 
perpendicular to local sweep and is positive outboard. 

Height in the boundary layer in DSTRZ units at which the 
crossflow velocity, in the plane perpendicular to the edge 
velocity, has its maximum absolute value. 

Height in the boundary layer in DSTRZ units at which the 
crossflow velocity is Xo of its maximum absolute value. When 
there are several points in the boundary layer where this is 
true, the largest height is printed here. 

Crossflow Reynolds number. The length unit is the height in 
boundary layer where the crossflow velocity is Xo of the maxi- 
mum, and the velocity unit is the maximum crossflow velocity. 

Boundary layer Reynolds number based on DSTRZ and edge 
velocity perpendicular to local sweep. 

Reynolds number used in the stability codes. This is local 
boundary' layer Reynolds number based on the local external 
velocity and height in the boundary layer to 99% of the exter- 
nal velocity. 

This is the magnitude of the local edge velocity in feet per 
second. 

Local static temperature at the edge of the boundary layer in 
degrees Rankine. 
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MACK STAG TEMP 
(TSTAG39 DEG K) 


Local stagnation temperature at the edge of the boundary 
layer in degrees Kelvin. 

LOCAL EDGE MACH Mach number of the local edge flow. 

ZME39) 

The next group of output contains the velocity and temperature profiles at this station and their 
first and second derivatives with respect to the local nondimensionalized “vertical" coordinate. 
ETA is the “vertical” distance nondimensionalized by the displacement thickness determined 
from the velocity profile in the vertical plane perpendicular to local sweep. U is the velocity com- 
ponent in this direction, and W is the component perpendicular to U. T is static temperature 
nondimensionalized by the edge static temperature. CROSS is the crossflow velocity determined 
in the plane perpendicular to the local edge velocity. 

The two groups of output described above are repeated for each station on the wing at which 
boundary layer information is transferred from the A411 code to the stability codes. 

The last set of output from LAMSD placed on the OP-file is nearly a repeat of what is on the 
BLGP-file that was explained earlier. However, the format is different and the CROSSM parame- 
ter that is placed onto the BLGP-file is omitted here. 

MKMOD3 Output on the OTSXQQZ File 

This file contains an expanded form of the output associated with the eigenvalue search for distur- 
bance wave angles less than 70 deg. The more brief form is on the OP-file. If a run fails in the 
MKMOD3 program, this output may be useful for finding and working around the problem. 

The first output on the OTS-file is a line with four numbers on it. This line is printed every time the 
program changes frequency in its solution process (i . e . , NFR times per station) . The first number 
is FR, the nondimensional frequency. The second number is PSI, the wave angle in degrees. The 
third number is the eigenvalue a, guess supplied to the solution process, and the fourth number is 
the Oj guess. 

The next set of output is repeated for each of the NPSI wave angles analyzed at the frequency 
printed on the preceding line. The words SPATIAL and COMPRESSIBLE remind the user of the 
type of solution process being used. In MKMOD3, only the SPATIAL COMPRESSIBLE equa- 
tions are allowed. The words NO DISSIPATION will appear next if the user has chosen to use the 
faster sixth-order equation set. ICNTRL and KMUL are two internal flags that the user need not 
be concerned about. ALPHA, DEL are the real and imaginary parts of the eigenvalue followed by 
a convergence parameter, and they are printed for each iteration in the solution process. The 

convergence parameter, DEL, in MKMOD3 is defined as / (Aa r ) + (Aa,) . Aoj r and Aoij are 

y a T 2 + at 

the predicted changes to 0t r and Oii for the next iteration in the Newton-Raphson search proce- 
dure; see section 4.2. The solution is considered converged in MKMOD3 when DEL is 0.005 or 
less. If the solution converges, the next set of output contains the following variables: 
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XQC 

ALPHARC 


ALPHAIC 


FR 


x/c at this station. 

a * r • c, this is the dimensional wave number in the unrotated 
coordinate system (see fig. 9) multiplied by the streamwise 
chord. 

a*i- c, dimensional spatial amplification rate in the unrotated 
coord, system (see fig. 9) multiplied by chord. This equals 
-(dN/ds) • c. 

2 

Nondimensional disturbance frequency,^* v* e /U* e . 


PSI3D 


Disturbance wave angle in degrees, . 


PSIB3D 


User-specified group velocity angle in degrees \j / . 


INTEGRAL 

BETARC, 

BETAIC 


UNROTATED 
ALPHA, BETA 

ALPHAR 


ALPHAI 


R 


Not used in MKMOD3. 

These are similar to ALPHARC and ALPHAIC but are in the 
“Z” direction in the unrotated coord, system. BETAIC will 
be 0. by definition if is input as 0. 

These are a r , a lt ft, and ft in the unrotated coord, system, ft 
is 0. by definition if \j? is input as 0. 

a r , nondimensional wave number in the rotated coord, system 
(see fig. 9). 

ai, nondimensional spatial amplification rate in the rotated 
coord, system. 

Re^, Reynolds number based on 8. 



Legend: 

X. Z Define the unrotated coordinate system. 

X, Z Define the rotated coordinate system, 
a Component of the wave number in the X or X 
direction. 

It has real and imaginary parts. 

(3 Component of the wave number in the Z or Z 
direction. 

It has real and imaginary parts. 


for Stability Analysis 

6-U90 150-0 
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CR 

Cl 

FR 

FI 

FREQ 

FREQI 

PHASE VEL 

PSI3D 

BETAR 

BETAI 


Real part of the phase velocity. 

Imaginary part of the phase velocity. 

Explained above. 

Imaginary part of complex frequency. This is 0. for this appli- 
cation of MKMOD3. 

FR times R. 

0. for this application. 

Dimensional phase velocity, FREQ/a r . 

Explained above. 

fi T in the rotated coordinate system, 0. by definition. 
fii in the rotated coordinate system. 


If the solution process does not converge, a message “MAXIMUM NUMBER OF ITERATIONS. 
CONVERGENCE NOT ACHIEVED” may be printed. If it is, the program will continue to run 
and may be able to successfully calculate eigenvalues at the conditions that follow. A more serious 
consequence of the solution not converging is an “arithmetic error failure, which will stop 
MKMOD3 calculations altogether. 


MKMOD3 Section of the OP-File 

A loadmap of the MKMOD3 program is the first output to be put on the OP-file having to do with 
the MKMOD3 program. This can be useful for locating serious errors that terminate the program 
during execution. 

The next group of output from MKMOD3, which appears on the OP-file, consists of two lines. It 
gives the wing station number (referenced to the stations picked by the LAMSD program), the 
station xJc location, the critical Reynolds number at this station (based on 8 ), the actual Re- 
ynolds number at this station (also based on 8*), and the Falkner-Skan p at this station. The 
velocity profile used to calculate 8* and the Falkner-Skan p is the one in the plane containing the 
velocity vector outside the boundary layer and perpendicular to the surface. 

If the program determines that this station may have disturbance growth, the next group of output 
gives the following information. The eigenvalue guess ( cir * &0 fur the first wave angle at this fre- 
quency is given first. This guess for the first eigenvalue is calculated by the main program and is 
determined by empirical knowledge and/or converged results at other frequencies or stations dur- 
ing this run. Guesses for the eigenvalues at the second and higher wave angles at this frequency are 
determined by logic in parts of the MKMOD3 code not modified under this contract. The next 
output is a r , a u and a convergence parameter called DEL and is printed for each iteration in the 
solution for this particular eigenvalue. If convergence is achieved, a message to that effect is 
printed next. If the solution does not converge, the program may terminate because of a numeri- 
cal error or a message will be printed that says “MAXIMUM NUMBER OF ITERATIONS, CON- 
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VERGENCE NOT ACHIEVED.” If this message appears, the program will at least try to continue 
the calculation process. Convergence is considered to take place when DEL is less than 0.005. 


When eigenvalues for each of the NPSI wave angles are calculated at the first frequency, a sum- 
mary is printed next which contains: 


R 

Re,, at this station. 

ALPHAR 

a r . 

ALPHAI 

a { . 

PHVEL 

This is zero for the mode being used in MKMOD3. 

FREQR 

Nondimensional frequency, FR times R. 

FR 

Nondimensional frequency. 

The program then increments the 

frequency and repeats the eigenvalue search through the wave 


angles again. 

After NFR frequencies are completed, the stability work at this station is complete and a summary 
table for this station is printed. The station number, x/c location, and angle of the local edge flow 
with respect to the freestream are printed first. The edge flow angle is in degrees and positive is a 
flow outboard with respect to the freestream. In the summary table the columns are at constant 
frequencies, which are in Hertz and are printed at the top of each column. The next row of 
printing gives the wave angle in degrees (referenced to the local edge flow direction, positive 
outboard) and then the disturbance growth rate, dN/ds, for this wave angle at each of the frequen- 
cies. The units for dN/ds are 1/ft. The third row of printing is the group velocity direction for the 
wave angle-frequency combinations. It is in degrees, referenced to the local edge flow direction 
and positive outboard. The above summary output continues until results for each of the NPSI 
wave angles are printed. The program then continues to the next station and the cycle of printing 
described above repeats until the run is complete. 

MKMOD3 Output on the TSIXQQZ File 

The TSI-file is used by the INTTS2 program to integrate the disturbance growth rates according to 
user-specified rules to get the N-factors at the various wing stations. The information on this file is 
a series of tables summarizing the MKMOD3 program results and is very similar to the tables that 
the program prints on the OP-file. 

The first line of output on the TSI-file contains, in order; NPSI, the number of different wave 
angles analyzed for each station; NFR, the number of different frequencies analyzed for each 
station; and the chord length in feet. The second line contains the x/c value of the station immedi- 
ately ahead of the first station at which the program starts calculations. If the program started 
stability calculations at the first station, then this value would be 0.0. 

The third line starts a set of output that is repeated for each station. This line contains, in order, 
the number of this station, the x/c value of this station, the s/c value of this station measured from 
the attachment line, the local sweep in degrees defined conventionally, and the angle in degrees 
between the local external velocity and a line perpendicular to the local sweep, defined as positive 
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for flow toward the wingtip. The second line of output in this set contains the values of the fre- 
quencies analyzed at this station. These are in Hertz and start with the lowest value at the left. 
After this line there are four lines, repeated for each wave angle. The first of these four lines starts 
with the disturbance wave angle in degrees and then has the disturbance growth rates, dN/ds, for 
that wave angle and the disturbance frequencies printed above it. The second of the four lines 
contains the nondimensional wave number, a,, which the program converged to for this wave 
angle and frequency. The third line contains the spanwise component of wave number in dimen- 
sional units of one per foot. Neither the dimensional nor non dimensional wave numbers in the 
TSI-file are presently used in the N-factor integration, but they are included for possible future 
use. The last of the four lines has the group velocity angle of this disturbance in degrees, measured 
relative to the local external velocity and positive outboard. After these four lines are repeated 
NPSI times, the summary information at this station is complete. The summary information for 
the next station is printed next and this continues for each station until the run is complete. 

MKMOD5 Output on the OCFXQQZ File 

This file contains an expanded form of the output associated with the eigenvalue search for distur- 
bance wave angles greater than 72 deg. The shorter form of this output is on the OP-file. If a run 
fails in the MKMOD5 program, this output may be useful for finding and working around the 
problem. 

The output on the OCF-file is similar to that on the OTS-file, which was explained earlier; there- 
fore, only the differences will be detailed here. For the MKMOD5 program the first line of output, 
which contains FR, PSI, the a r guess, and the a, guess, is printed for each eigenvalue solution, 
instead of only once for each different frequency, as was the case with the MKMOD3 program. 
After the ICNTRL, KMUL line there is a message printed: VARYING BETAI — NOT ALPHAI. 
This tells the user that the program is operating at high-disturbance wave angles, greater than 72 
degrees. Before the ALPHA, DEL line of output, the OCF-file has a line containing DELT1, 
DELT2, and AMPAZT. DELT1 is the change in nondimensional frequency predicted by the 
program for the next iteration. DELT2 is the change predicted for f $\ , which is the nondimension- 
al dN/ds for MKMOD5. It is important to note here that the eigenvalues solved for by MKMOD5 
during the iteration process are frequency and fa , with wave angle, \|r. and wave number, a r , held 
fixed. In the MKMOD3 program, a T and are the eigenvalues solved for, with frequency and 
wave angle held fixed. The AMPAZT parameter gives the user some idea of how well the solution 
process is matching the boundary conditions at the wing surface. The reader will recall that the 
basis for this stability method is an integration of the disturbance equations through the boundary 
layer from top to bottom, changing the eigenvalues until the wall boundary conditions are ade- 
quately met. The magnitude of AMPAZT is not of much concern to the user, but it should de- 
crease as the iteration process converges. When the solution is converged, AMPAZT will 
generally be 10 or less, but for some very difficult velocity profiles, AMPAZT may be a very strong 
function of frequency and fa . In these cases AMPAZT may still be in the order of 1000 when 
sufficient eigenvalue accuracy has been achieved. 

Following the line described above is the ALPHA, DEL line. As mentioned before, the real part 
of ALPHA is held fixed in this program, although the imaginary part varies. In MKMOD5, DEL is 
defined as 


Jl . \ 2 , 

( DELT2 In dimensional units \ * 

V \ 1 X 10-6 J 

\ 5. ) 
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For convergence DEL must be less than or equal to 0.01. When convergence occurs, the output 

that starts with the “XQC, ALPHARC ” line and ends with the “CONVERGENCE 

ACHIEVED" line is printed and is the same as was explained for the MKMOD3 program. 

The output explained above is repeated for each eigenvalue that the MKMOD5 program finds or 
attempts to find. 

MKMOD5 Output on the OPXQQYZ File 

As with the MKMOD3 program, a loadmap of MKMOD5 is the first output that is printed on the 
OP-file. 

The next set of output from the MKMOD5 program that appears on the OP-file has two lines. 
Station number and x/c are given on the first line. The station refers here to the stations picked by 
the LAMSD program. The second line contains the critical crossflow Reynolds number, the actual 
crossflow Reynolds number at this station, and the value of the crossflow H parameter. The criti- 
cal crossflow Reynolds number is a function of crossflow H and is determined from an empirical 
table lookup from reference 30 (see fig. 8). Only if the crossflow Reynolds number is 10 or more 
units greater than the critical crossflow Reynolds number will the eigenvalue calculations begin. 

If eigenvalue calculations start at this station, the next output is the eigenvalue guesses at the first 
wave number/wave angle calculation point. The wave number is chosen from expirical informa- 
tion by the program. The wave angle is the so-called “critical angle” minus a user-supplied incre- 
ment, DPSIST, generally between 1 and 3 deg. The ALPHAR printed right after EIGENVALUE 
GUESS is the fixed wave number to be used at this condition and the BETAI and FR are the 
guesses for the two eigenvalues. At the first conditions these are the user-supplied inputs called 
ALPHAI and FRINP, respectively (see sec. 6.0) . Results of the iteration process are printed next. 
The first number to the right of ALPHA, DEL is a,. This remains fixed in the MKMOD5 iteration 
process. The next number is a,, and it changes as #is changed during the iteration. DEL is the 
measure of the change in the eigenvalues requested by the program for the next iteration. When 
DEL is 0.01 or less, the solution is considered converged. If the solution converges, a message to 
that effect is printed as well as the final values of the eigenvalues. The program then goes to the 
next wave angle/wave number condition, calculates the eigenvalue guesses for that condition, and 
repeats the solution process. For this first series of solutions the program is finding the lower limit 
on the wave angle range to use for this first station. It does this by fixing the wave number and 
finding solutions at different wave angles, starting at the first angle and working down in 0.5-deg 
increments. This series continues until either the wave angle is less than 73 deg or a converged 
BETAI is found that is less than 30% of the first BETAI. If the eigenvalue iteration process does 
not converge, the eigenvalues that are returned from it will be much in error and the wave angle 
ranging logic will probably fail. An improved user-input guess may solve this problem. If the wave 
angle range is calculated normally, this series of solutions will be terminated with a printout saying, 
WAVE ANG. RANGE FOR FIRST STA. IS _ TO The maximum value for the 

wave angle is always the “critical angle” plus 0.6 deg. 

Once the wave angle range for the first station is determined, the wave number range is found. 
This is done by going back to the first wave angle and decreasing the wave number by 0. 1 for a 
series of conditions that end when 10 different wave numbers have been analyzed at this wave 
angle or until a converged BETAI is found that is less than 20% of the first BETAI found. The 
program then switches from decreasing the wave number to increasing it from the very first condi- 
tion. The results of this series of wave number calculations are curve fit, and the wave numbers to 
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analyze for the first station are determined from it. A message is printed giving the wave number 
range. The range printed here applies to the (NPSI-1) wave angle only, because the wave number 
range changes with wave angle. Experience with MKMOD5 has shown that the unstable region of 
wave angle/wave number conditions can best be covered for a wide variety of boundary layer 
profiles if the lower limit for wave number is held fixed for the different wave angles, but the upper 
limit is a function of wave angle to the 1.6 power. 


After the ranges of wave angles and wave numbers to analyze are found for the first station, the 
standard order of calculating eigenvalues for these ranges is worked through. The program starts 
at the lowest wave angle and wave number and proceeds from there to higher wave numbers. 
When all the NFR wave numbers are analyzed at the first wave angle, the calculation goes to the 
next higher wave angle and starts again at the lowest wave number, working to higher ones. The 
output seen by the user on the OP-file during this process is a set of EIGENVALUE GUESS — , 
ALPHA, DEL, and RESULTS FROM MACK, for each converged solution. If the solution does 
not converge during this part of the analyses, it will trigger internal recovery logic which will either 
change the eigenvalue guess (this is the approach used for the first four wave angle and wave 
number combinations) or break up the step from a previous converged solution into 8 small steps 
so the eigenvalue guess can be improved. This recovery procedure is explained more fully in figure 
10 . 
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Figure 10. — Estimating the Next Eigenvalue and Recovering From Divergence — MKMOD5 
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When eigenvalues for all the wave angle/wave number combinations for this station have been 
calculated, a summary of the results is printed. The summary contains the station number and x/c, 
the critical angle in degrees as well as the angle of the local edge velocity with respect to the 
freestream, also in degrees and positive outboard. The next five lines give calculation results at the 
lowest wave angle. Included here are wave angle in degrees, frequencies in Hertz, nondimensional 
wave numbers, the dimensional wave number components perpendicular to the local edge veloc- 
ity; the disturbance growth rates per foot, dN/ds; and the group velocity angles of these distur- 
bances in degrees, referenced to the local edge velocity and positive outboard. This set of five 
lines is repeated for each of the wave angles analyzed. 

The output for this station is completed with the printing of the wave angle ranges to be used for 
the next station. These are calculated by curve-fitting with parabolas the dN/ds versus wave num- 
ber, a r , at the first and next-to-highest (NPSI-1) wave angles. If the maximum value of dN/ds of 
this parabola is greater than 1.0, the wave number range for the next station is determined to be 
from the wave number where dN/ds — 0 on the lower branch of the parabola to the wave number 
where dN/ds a negative % of the maximum value on the upper branch. If the maximum value of 
dN/ds of the parabola is less than 1.0, the wave number range is from the wave number on the 
lower branch at which dN/ds equals the maximum value of dN/ds on the parabola minus 1 .0 to the 
wave number on the upper branch at which dN/ds is the maximum value minus 1.5. In addition, 
the lower end of the range of wave number must be greater than or equal to 0.3. In this way the 
wave number ranges for all stations after the first one analyzed are determined by the eigenvalue 
results from the previous station. 

As with the wave number range determination, the lower limit on the wave angle range for all 
stations after the first is calculated using results from the previous station. Again a curve-fit is used. 
In this case dN/ds versus wave angle results for the (NFR+l)/2 wave numbers are used. In the 
summary table described earlier, this is the middle column for odd NFR or the column just to the 
left of center for an NFR that is even. The lower range of wave angle is then determined to be at 
the angle where dN/ds is % of the maximum value from the parabolic curve-fit, or 72 deg, which- 
ever is larger. The largest wave angle to be considered at the next station is again the critical angle 
for that station plus 0.6 deg. The wave angle range for all stations after the first one analyzed is 
printed after the first two lines of output for that station. 

The Qutput described above that did not have to do with finding the wave angle and number 
ranges for the first station analyzed is repeated in sets for the rest of a MKMOD5 run. 

MKMOD5 Output on the CFIXQQZ File 

The CFI-file is similar to the TSI-file in that it summarizes the MKMOD5 program results and is 
used by an integration program, INTCF2, to calculate N-factors. The series of tables printed on 
the CFI-file are similar to those that MKMOD5 prints on the OP-file. 

The first line of output on the CFI-file contains, in order; NPSI, the number of different wave 
angles analyzed during this run; NFR, the number of different wave numbers analyzed during this 
run; and the chord length in feet. The second line contains the x/c value of the station immediate- 
ly ahead of the first station at which the program starts calculation. If the program started stability 
calculations at the first station, then this value would be equal to x/c), - [x/c) 2 - x/c),] or 0., 
whichever is greater. 
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The third line starts a set of output that is repeated for each station. This line contains, in order, 
the number of the station, the x/c value of this station, the s/c value of this station measured from 
the attachment line, the local sweep in degrees defined conventionally, and the angle in degrees 
between the local external velocity and a line perpendicular to the local sweep, defined as positive 
for flow toward the wingtip. The fourth card is the beginning of a set of five cards that summarize 
the MKMOD5 results for the first wave angle examined at this station. This line contains first the 
wave angle in degrees and then all the disturbance frequencies for each of the wave numbers 
examined at this wave angle. The frequencies are in Hertz and correspond with the nondimen- 
sional wave numbers that are printed on the line immediately following them. The third line in this 
set is the component of the wave numbers perpendicular to the local edge velocity in dimensions 
of one per foot. The fourth line contains the disturbance growth rates in dimensions of one per 
foot for each wave number printed immediately above. The fifth line gives the corresponding 
group velocity angles in degrees, measured with respect to the local edge velocity and positive 
outboard. The output then proceeds to the next wave angle and repeats the above set for it, 
continuing through all the wave angles examined at this station. The same information is then 
printed for following stations until the run is complete. 

INTTS2 Output on the OTSIXQQ File 

The INTTS2 program creates two output files: one is readable and the other is for machine plot- 
ting. Both files contain N-factors as a function of x/c for different frequencies. The INTTS2 
program presently calculates N-factors by both the maximum amplification method and the meth- 
od used by Boeing (constant wave angle method) . Results from both these methods are printed on 
the OTSI-file, starting with the maximum amplification method. The frequency in Hertz for the 
first N-factor integration is printed first. Following that, four columns are printed that are headed 
with their descriptions: x/c, dN/ds, wave angle, and N-factor. The program searches the table of 
dN/ds versus frequency and wave angle for each station (TSI-file) to find the wave angle for 
maximum growth rate at the desired frequency. This procedure is allowed to extrapolate outside 
the table but the resulting value of dN/ds found from extrapolation may not be reasonable. There- 
fore, if extrapolation occurs, a message to that effect is printed and the user should closely ex- 
amine the associated N-factors. The next set of output is for the next frequency analyzed. Note 
that the program searches all the tables contained in the TSI-file to find the frequency range over 
which to calculate the N-factors, and then goes back and calculates N-factors for 36 different 
frequencies in this range. This range of frequencies has always proven to be more than adequate 
to define the envelope of most amplified disturbances. 

After N-factor results for the 36 frequencies using the maximum amplification method are 
printed, the program repeats the process for the same frequencies, using the first of the user-speci- 
fied wave angles. This is the method used by Boeing. Several different wave angles can be speci- 
fied, and N-factors for each are printed as a set in the order in which they are input (see sec. 6.0) . 

INTTS2 Output to the GTSIXQQ File 

The GTSI-file contains N-factor information from the 1NTTS2 program in plotting format. The 
first set of output before an end-of-file mark are the results from the N-factor calculations using 
the maximum amplification method. The first line contains the format of each dataset. The sec- 
ond line is a graph title. The third line is a note placed on the graph which tells the reader that the 
NTS - x/c plot was calculated using dN/ds at the wave angle that gave maximum growth. Lines 
four through seven give first the independent variable name, XOC, and then the parameter names 
associated with each of the N-factor versus x/c plots that can be made. These parameter names 
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are “F” followed by the frequency in Hertz rounded to the nearest whole number. There is one 
parameter name for each of the 36 frequencies at which the INTTS2 program calculated N-fac- 
tors. What follows this line of output are sets of numbers, one for each station at which MKMOD3 
analyzed the boundary layer stability, starting at the first station. The first number in each of these 
sets is the x/c of the station. The numbers that follow are the N-factor information for the NASA 
calculation method. 

For N-factor calculations using the constant wave angle method, a set of output almost identical to 
that described above is created for each wave angle specified by the user. The graph title changes 
appropriately, as does the “note” given on line three. Each set ends with an end-of-file mark. 

INTCF2 Output on the OCFIXQQ File 

The INTCF2 program creates two output files: one for reading and the other for machine plotting. 
Both contain N-factors as a function of the x/c for (1) different frequencies, in the case of the 
maximum amplification method, (2) different values of dimensional components of wave number 
in the spanwise (perpendicular to local edge flow) direction, in the case of the irrotational integra- 
tion method. Results from both these methods are printed on the OCFI-file, starting with the 
maximum amplification method. The frequency in Hertz for the first N-factor integration in 
primed first. Following that, four columns are printed that are headed with their descriptions; x/c, 
dN/ds, wave angle, and N-factor. The program searches the table of dN/ds versus frequency 
(which varies with wave number) and wave angle for each station to find the wave angle for 
maximum growth rate at the desired frequency. If the program determines that the wave angle for 
maximum growth rate at this frequency is at a wave angle that is lower than that available in the 
CFI-file, it includes information from the lower angle analyses on the TSI-file, providing that file is 
in existence and has information on it at the desired station. The maximum amplification integra- 
tion method is allowed to extrapolate beyond its tables in some limited ways. If it does extrapolate 
or pick growth rate values at the edge of the tables, it prints a message that warns the user that the 
associated N-factors may not be accurate. The next set of output is for the next frequency. The 
program searches all the tables contained in the CFI-file, but not the TSI-file, to determine the 
frequency range over which to calculate the N-factors. It then calculates N-factors for 25 different 
frequencies in this range, which usually is adequate to define the envelope of most amplified 
disturbances. Since this integration program can use disturbance growth rates from wave angles 
calculated by both of the stability programs, the N-factors calculated by the maximum amplifica- 
tion method from this program are probably more accurate than those calculated by the INTTS2 
program. 

After N-factor results for the maximum amplification method are printed, the program switches to 
the irrotational method. After the printed message to this effect, the frequency in Hertz and 
spanwise wave number in units of 1 per foot are printed. Following this, the four columns of x/c, 
dN/ds, wave angle, and NCF are printed for this wave number. Extrapolation beyond the wave 
angle and frequency tables provided by the CFI-file is also permitted for this integration method, 
and a message warning the user is printed if extrapolation is done. This set of output is repeated 
for each of 26 different wave numbers chosen automatically by the program. If the user has 
specified by input to the INTCF2 that N-factors for more than one frequency are to be calculated 
with the irrotational method, the above output is repeated in turn for each frequency. 
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INTCF2 Output to the GCFIXQQ File 


The GCFI-file contains N-factor information from the INTCF2 program in plotting format. The 
first set of output before an end-of-file mark are the results from the N-factor calculations using 
the maximum amplification method. The first line contains the format of each dataset. The sec- 
ond line is a graph title. The third line is a note placed on the graph that tells the reader the 
N - x/c plot was calculated using dN/ds at the wave angle that gave maximum growth. Lines four 
through six give first the independent variable name, XOC, and then the parameter names asso- 
ciated with each of the N-factor versus x/c plots that can be made. These parameter names are 
“F” followed by the frequency in Hertz rounded to the nearest whole number. There is one 
parameter name for each of the 25 frequencies at which the INTCF2 program calculated N-fac- 
tors. What follows this line of output are sets of numbers, one for each station at which MKMOD5 
analyzed the boundary layer stability, starting at the first station. The first number in each of these 
sets is the x/c of the station. The numbers that follow are the N-factors at this station for each of 
the frequencies, and in the same order as the parameter names on lines 4 through 6. The set of 
N-factor information for the maximum amplification calculation method ends with an end-of-file 
mark. 

For crossflow N-factor calculations using the irrotational method a set of output very similar to 
that described above is created for each frequency specified by the user. The graph title changes 
appropriately, as does the “note” given on line three. The parameter names are “K” followed by 
the spanwise component of wave number rounded to the nearest whole number. The set of N-fac- 
tor information for each frequency is separated by an end-of-file mark. 
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8.0 GUIDELINES FOR THE USS USER 


8.1 PRERUN CONSIDERATIONS 

Preparing to make a run on the Unified Stability System (USS) starts with gathering the informa- 
tion needed in the input deck. Of the information required, the accuracy of the C p - x/c distribu- 
tion is among the most critical. The more sparse the definition of the pressure distribution or 
questionable its accuracy, the more care will be required to ensure that the dense definition gener- 
ated by the BLGL program for the boundary layer analysis is a “reasonable” one. Reasonable as 
used here means that aerodynamic judgment may be required. The user may be disappointed to 
see that the pressure distribution for the complete wing section is required as input. Although this 
is the case, the accuracy of the geometry and pressure outside the region of interest is not very 
important, except as it may affect the curve-fitting near the edges of the region of interest. This 
may be an important issue near the leading edge when pressure data are only available on one 
surface or the other. 

Another important aerodynamic input is the taper (the actual input is local sweep at two different 
chordwise locations) . If nothing is known about the isobar pattern in the area of interest, geomet- 
ric taper can be used. However, the boundary layer really responds to aerodynamic taper, so if the 
aerodynamic taper changes with x/c in the region of interest, some average must be chosen. As a 
rough rule of thumb, the isobars that the USS assumes based on its two sweep inputs, should 
match the actual isobars most closely in regions of high chordwise pressure gradient. 

8.2 CONSIDERATIONS FOR BOUNDARY LAYER CALCULATIONS 

When a new USS run is started, the user can choose to run just the first program (BLGL) or as 
many as all seven of the programs in one job submittal. If there is a question concerning the 
curve-fitting of the pressure distribution by BLGL, only that program should be run and the dense 
array of C p - x/c that it generates should be plotted. A distribution that results in a poor curve fit 
can be improved by adding additional input points. If the curve-fit of pressure seems adequate, 
the second program in the USS, the A411 boundary layer code, can be run. It is possible that 
A41 1 will predict laminar separation and stop calculations at a point on the wing section ahead of 
where the user knows transition took place. If this happens, all input that affects the boundary 
layer calculation should be checked first. If the input is correct and the pressure distribution 
curve-fit looked acceptable, the user has little recourse but to make small changes to the pressure 
coefficient input until the program runs at least to the point where the user knows transition took 
place. 

After the boundary layer is calculated, appropriate stations need to be picked for input into the 
stability programs. This is done by the third code of the USS, LAMSD. As described in the input 
description, there are three options for choosing the stations for stability analyses from the many 
calculated by A4 1 1 . Of these, the option in which the user specifies the s/c locations of the desired 
stations is recommended. The distribution and number of stations to be picked varies consider- 
ably with each case, but there are guidelines that can assist the choices. In general, pick stations 
about 0.5% chord apart for the first 5% chord, starting at 0.5% chord. For the next 5% to 10% 
chord, pick stations about 1% chord apart. For the rest of the distance, 2.5% to 5.0% chord 
spacing may be possible. Without knowing the boundary layer results, a priori, it is the pressure 
gradient and its variation that guides the user in determining the station distribution. Rapid 
changes in gradient call for closer spacing. Closer spacing also implies smaller changes in stability 
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characteristics between stations, which makes the eigenvalue “guessing” logic in the stability pro- 
grams more accurate. The number of stations chosen is a compromise between wanting a detailed 
picture of the stability characteristics and the cost of the computation time. There is also a limit of 
40 stations currently in the USS. 

The LAMSD program can also create a file containing boundary layer profiles in a plottable 
format. These are generally not required, but some cases may exhibit such strange stability behav- 
ior that the user will want to go back and rerun LAMSD to see the profiles. The means of picking 
the stations for which profiles will be saved on the plotting file is cumbersome (see the input 
description) , but does allow the user to plot the profiles of any station calculated by the boundary 
layer program. The correlation between A411 station number, x/c at the station, and s/c at the 
station can be found on the CREF-file created by the BLGL program. 

8.3 SETTING UP THE MKMOD3 RUN 

The MKMOD3 program calculates the boundary layer stability characteristics at wave angles be- 
low 65 deg. It develops sets of information that can be represented as tables as shown in table 1. 


Table 1. — MKMOD3 Program Tables 


Wave angle 

Frequency, Hz 


• Wj 

W NFR 


dN > 

ds J 

h. 


Disturbance growth rate 
Group velocity angle, deg 

dN\ 

dS / NFR, 1 
'/'NFR, 1 

V'K 







'/'NPSI 





dN /ds) NFR NPS , 
'/'NFR, NPSI 


For the first stability analysis run of a given case, the user usually specifies that the program start at 
station 1 and continue to the last station. The program then determines from the boundary layer 
characteristics at which station there may be unstable disturbances and starts eigenvalue calcula- 
tions there. This logic is based on comparing the displacement thickness Reynolds number of the 
boundary layer profile perpendicular to the local sweep with a “critical” Reynolds number from 
reference 29 (see fig. 4) . Once the program is triggered to start eigenvalue calculations, it does not 
“turn off” at later stations that may have no disturbance growth. 

The DELPSI input to MKMOD3 determines how much of the wave angle region is covered by this 
run. Internal limits constrain the minimum wave angle to be -50 deg and the maximum to be 65 
deg. Picking DELPSI is a compromise between the desire to— 

a. Establish the stability characteristics over a wide wave angle range. This is particularly impor- 
tant for using the “maximum amplification” N-factor calculation. 

b. Limit cost by choosing fewer wave angles to analyze, and/or the desire to get a more dense 
definition of the stability with a given number of wave angles, and/or the desire to increase the 
odds for a successful run by making the increments between wave angles smaller. 
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The program tries to distribute the wave angles it analyzes around an empirical estimate of the 
most unstable wave angle, which is a function of local Mach number (see ref. 6). In practice this 
empirical estimate of the most unstable wave angle has not proven to be very accurate, but it is still 
worth including. As an example, if DELPSI = 20 and NPSI = 5, the wave angles of about -40, -20, 
0, 20, and 40 would be analyzed for local Mach numbers below 0.7. For local Mach numbers in 
the 1.0 to 1.3 range, the upper limit on wave angle would come into play. The lowest angle would 
be near 0. deg. The resulting wave angles to be analyzed would then be approximately 0, 16, 33, 
49, and 65 deg. 

The frequency range to be analyzed by MKMOD3 is chosen by empirical logic in the program, 
developed from reference 29 (see sec. 4.2). For most cases this range adequately covers the 
region of instability, at least well enough to establish the envelope of disturbance growth. The 
empiricism is weakest for boundary layers that are near separation. 

8.4 CHECKING FOR PROBLEMS AND RESTARTS WITH MKMOD3 

The tables of disturbance growth generated by MKMOD3 are available on the OP-file or the 
TSI-file. It is advisable for the user to scan these tables to check for unusual results. The program 
may occasionally return incorrect results, even if the failure is not serious enough to stop the run. 
One possible failure is nonconvergence of an eigenvalue search that is not serious enough to stop 
the run. This can be spotted in the stability tables if dN/ds at two or more frequencies have 
identical values for the same wave angle or if the group velocity angles appear to be too high. The 
group velocity angles are calculated by an interpolation process. This is an efficient but not a 
highly accurate method, and the answers from it are sometimes suspicious; occasionally they are 
obviously wrong. This is not necessarily a cause for concern. Although the group velocity angle is 
used in the N-factor integrations, the integration programs limit the angle they use to be between 
±10 deg. If the stability tables have strange-looking results, the user should look at the OP-file 
and perhaps the OTS-file to get more clues as to the validity of the results. 

If the user wishes to start the MKMOD3 program at some station partly through the region of 
interest, this presents no particular problem. The start-up logic for guessing eigenvalues at the first 
station calculated in a run is sometimes even more robust than that which uses results at previous 
stations. This is an especially important point if the boundary layer changes are excessive between 
stations. On a run that continues calculations begun in a previous run, the user must change the 
version number or the OP-, OTS-, and TSI-files created by the second run will write over those 
created by the previous run. For a case that has been completed by several different runs, the user 
must combine the TSI-files from the runs in order for the integration programs to work properly. 
Additionally, NFR and NPSI for all these files must not change. 

If the MKMOD3 program fails, consider the following: 

a. Density of the stations with regard to the change in the boundary layer between adjacent 
stations. 

b. Size of the wave angle increment. 

c. Size of the frequency increment. 
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8.5 SETTING UP AND RESTARTING AN MKM0D5 RUN 

The MKMOD5 program calculates the higher wave angle stability characteristics. It has a lower 
limit on wave angle of 72 deg and an upper limit of 0.6 deg above the critical angle of the 
boundary layer profile. This wave angle range does not necessarily cover all the unstable area but 
is intended to cover all the unstable area for positive frequencies. For most cases of interest to 
wing designers, the envelope of the most amplified disturbances will be defined by positive fre- 
quency disturbances. 

Like the MKMOD3 program the MKMOD5 program develops “tables” of stability information at 
each wing station. 


Table 2. — MKMOD5 Program Tables 
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The frequency in this program is not constant on a column of the table, as it was for the MKMOD3 
program. 

If the MKMOD5 program is started from station 1 for a stability analysis run, it parallels the 
MKMOD3 program in examining the boundary layer characteristics to determine which station 
probably contains disturbance growth. In MKMOD5 this logic compares the crossflow Reynolds 
number of the boundary layer profile to an empirically determined critical value of crossflow 
Reynolds number as given in reference 30. This program, unlike MKMOD3, does contain logic to 
skip stability analyses at stations if (1) the maximum dN/ds at the (NPSI-1) wave angle of the 
previous station was less than -2. and (2) the crossflow Reynolds number is 15 less than the 
empirical “critical” value. At the high wave angles the solution for the eigenvalues is more difficult 
than at lower wave angles, especially for conditions in which the disturbances are damped. This is 
why logic in the MKMOD5 program is tailored to prevent calculations at stations that are probably 
stable. A compromise on this issue is required; however, because an adequate definition of distur- 
bance growth and damping is required to properly define the N-factor envelopes. 

For cases showing a “reasonable” amount of crossflow instability near the leading edge, the sug- 
gested input values for ALPHAI, FRINP, and DPSIST will probably result in a normal program 
start. Cases that are only slightly unstable but do trigger the logic that starts eigenvalue calculations 
may require changes in ALPHAI and FRINP to obtain program convergence during the starting 
phase. This may also be true when starting the program at downstream stations. The starting phase 
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involves finding the wave angle and wave number ranges to cover in subsequent table calcula- 
tions. After the first station these ranges are determined from results at the previous, upstream 
station. If the program does not get good converged eigenvalues during this starting phase, it does 
not necessarily fail, but the subsequent stability table results should be closely examined. The 
degree of instability at the starting station can be critical for getting a successful start; generally the 
more unstable the more likelihood of a successful start. Some measure of the instability at a 
station can be gained by comparing the crossflow Reynolds number to the critical crossflow 
Reynolds number. The former is located on LAMSD output to the OP- or BLGP-files, and the 
latter can be found from figure 8, knowing H (crossflow) . H(crossflow) is also output by LAMSD. 
The user’s desire to start the calculations at a station that is unstable enough to get a successful 
start must be tempered by remembering that N-factor accuracy is reduced when the chordwise 
increment between stations is increased. 

The stability tables generated by MKMOD5 may be found on the OP- or CFI-files. It is advisable 
to scan these tables to check for unusual results. After a few successful runs users will get an idea 
of what constitutes an “unusual” result. The program can give incorrect results without actually 
stopping. One possible problem is the inability to converge during an eigenvalue solution. Because 
this is much more likely to happen at these high wave angles, MKMOD5 has logic for recovering 
from cases that do not converge. For the first four points in the stability table (two lowest wave 
angles and wave numbers) , this logic involves changing the ALPHAI guess and trying the solution 
again, for a maximum of ten attempts. For the rest of the points in the table, the recovery logic 
involves taking a smaller step (% the size) in wave angle or wave number from the previous point, 
at which it is assumed there is a converged solution. Figure 10 outlines this recovery procedure. 
The group velocity angles are calculated by MKMOD5 the same way as done by MKMOD3, so 
the shortcomings mentioned with regard to that program are also possible for this one. If the 
stability tables have strange-looking results, the user should look at the OP-file and perhaps the 
OCF-file to get clues as to the validity of the results and possible causes of problems. 

Restarts of MKMOD5 are more common than with MKMOD3 because unusual program termina- 
tions are more common for it. Unfortunately, restarts are also more difficult with MKMOD5 
because— 

a. The logic affecting the start of calculations is less robust in MKMOD5 because the calculations 
are inherently more difficult at the high-wave angles. 

b. The logic that allows the program to start calculations at a station based on crossflow Reynolds 
number may not allow a restart in an area that is only slightly unstable. 

c. The initial eigenvalue guesses that the user must input may be more difficult to determine. If 
results at previous stations are available, they are usually helpful. 

8.6 DISTURBANCE GROWTH CALCULATIONS 

The two programs that use the stability tables to calculate N-factors are similar and will be dis- 
cussed together. Both programs are presently set up to do two sets of calculations, one using the 
constant wave angle (for TS) or irrotational (for CF) approach and one using the maximum 
amplification approach. The maximum amplification approach may require stability information 
over a wide wave-angle range, so if both TSI- and CFI-files are available, the INTCF2 program 
uses them both. However, the other integration program, INTTS2, still does the maximum ampli- 
fication integration as well because it is possible that the CFI-file might not have been calculated. 
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The user is reminded that if both files are available, the maximum amplification results from 
INTCF2 are probably more accurate. The following points should be remembered about the 
integration programs: 


a. Accurate N-factor calculations require knowing the disturbance growth rates from the point 
on the wing where they first become unstable. 

b. If the TSI- and CFI-files were generated by combining the results of several different runs of 
the stability codes, be sure there are no duplicate stations in the files. 

c. A failure of the integration programs is unusual. If it happens, check the CFI- or TSI-file and 
the tape 5 input. If there is a suspicious station in the CFI- or TSI-file, it may be edited from 
the file if the user feels the integration without that station will still have enough accuracy. 

d. The message printed by the integration programs concerning extrapolation is only informative . 
Extrapolation is permitted because in many cases, it improves the answer. However, whether 
it is reasonable or not in a particular instance has to be judged by the user. 

e. The boundary layer and stability characteristics are calculated on a streamwise section of a 
wing. However, the integration of dN/ds over the wing does include the effect on A s of the 
disturbances growing in the direction of the group velocity (see fig. 11). 
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Figure 11. — Effect on As of the Disturbances Growing in the Direction of the Group Velocity 
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9.0 CONCLUSIONS AND RECOMMENDATIONS 


The Unified Stability System (USS) was developed to provide a more versatile sweptwing laminar 
boundary layer stability analysis method and also one that was easier to use than previous meth- 
ods. These two goals tend to be at odds with each other, and obviously, compromises had to be 
made. Most of the computer coding necessary for the system, the laminar boundary layer analysis 
and the stability analysis, was a modification of existing codes, which saved a great deal of effort. 
Additionally, results from previous authors were very helpful in setting up the initial eigenvalue 
"guessing” logic used in the stability codes of the USS. 

The USS was thoroughly verified during its development. Some of this verification is described in 
reference 4 and some was a result of the F-l 1 1 reanalyses described in the appendix of this report. 

The numerous files and large amount of output generated by the USS can be intimidating for the 
first-time user. Most of this information is never used, but it is still included for possible debugging 
of problems or for use by the analyst looking at details of the boundary layer or its stability charac- 
teristics. Therefore, the cost of running a case (one streamwise “cut” on a wing) will depend on 
the computer system and costing algorithm used, but the USS can be considered a moderately 
expensive system to use. The first-time user of the USS should either carefully read the sections of 
this report on the system’s use or use the system under the direction of someone who is familiar 
with it. 

Several files created by a USS run are information formatted for plotting on Boeing graphics 
facilities. It is hoped that these files will require little change to make them available for plotting 
by other computer systems. 

Currently, the USS can be used for examining the stability characteristics of laminar boundary 
layers on sweptwings from incompressible speeds through transonic speeds. The extension of the 
USS to supersonic speeds is believed possible, but has not been verified at this time. 


The N-factors at transition presented in the appendix for F-14 VSTFE, F-l 1 1, and 757 NLF flight 
data show considerable scatter. This indicates a need for further development of correlations of 
transition on sweptwings. The USS should be a valuable tool for this research, eliminating much of 
the handwork previously required. 

It was shown in references 36 and 37 that the curvature of the path of the disturbance can intro- 
duce terms into the classical stability equations that noticeably affect the eigenvalues of these 
equations for the spinning disk and yawed cylinder cases respectively. Recently, the effect of these 
curvature terms has been examined for sweptwing boundary layers and has been found to have a 
substantial effect on stability there too, reference 38. It would be straightforward to add curvature 
terms arising from surface curvature and streamline curvature to the USS, but considering the 
curvature of the path of the disturbance followed by any particular growth philosophy would re- 
quire an iterative process. 
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11.0 APPENDIX-STABILITY ANALYSES USING THE USS 


Fifteen cases were selected from the Variable Sweep Transition Flight Experiment for boundary 
layer stability analysis. The measured data were obtained from the F-14 cleanup glove. The data 
for each case consisted of flight conditions, pressure distributions and hot-film traces on the upper 
surface. The results from these analyses provide further insight into how the crossflow (CF) and 
Tollmien-Schlichting (TS) disturbances interact during the boundary layer transition process, for 
a range of sweep angles. Table A-l lists the features of the 15 cases. 

Three separate rows of pressure taps were placed in the cleanup glove to measure the pressure 
distributions at wing-buttock-lines 200.8, 260, and 320 (defined with Ale at 20 deg). Each row 
consisted of 19 ports on the upper surface and 2 ports on the lower surface. In addition, the state 
of the boundary layer was determined by surface hot-film sensors. These sensors were installed at 
every 10% chord location, staggered about 30 deg from the streamwise line. Only one row of 
hot-film sensors was used per case, placed between the first and the second or the second and 
third pressure orifice rows. Figure A1 shows a schematic planform of the test region. 

For the selected cases, the three rows of pressure data were used to determine the pressure iso- 
bars. Knowing the leading edge sweep and isobars, an effective taper ratio for the laminar portion 
of the wing could be determined. The USS code has the capability to include the taper effect. The 
hot-film traces were used to locate the transitional region by locating the gauge that shows transi- 
tional traces or had just gone turbulent. The location of that gauge was then marked and an 
effective chord line was drawn through it. Figure A2 illustrates the definition of effective chord as 
opposed to chord at a nominal sweep angle. The pressures along the effective chord line w'ere 
found by an interpolation between the three measured rows to the WBL of the measured transi- 
tion point. This chordwise distribution was used for the boundary layer stability analysis. 


Once all the data were reduced to the final form required by the USS, the code was run for all 15 
cases. The compressible CF and TS disturbance growth envelopes, with the corresponding pres- 
sure distributions, are shown on figures A3 through A22. For TS disturbances, the solid line 
represents the envelope of growth of disturbances of different frequencies at a wave angle that 
gives the maximum growth. This is the “constant wave angle” approach described in the body of 
this report. However, if the growth envelope keeps increasing with wave angle, a limit of 55 deg is 
imposed. For many laminar boundary layers, there is no clear division between a region of TS 
disturbances and CF disturbances. For CF disturbances, the dashed line represents the envelope 
of crossflow disturbance growth at zero frequency and with a range of values of the spanwise 
component of the dimensional wavenumber, (this is the “irrotational approach). For a few 

of these cases, the “maximum amplification” envelope is also shown in an additional figure. That 
envelope represents the maximum growth of disturbances of varying frequencies, with wave angle 
changing to give maximum growth rate for each frequency at each chordwise position as the 
disturbance moves aft on the wing. The transition location for each case, determined from hot- 
film traces, was used along with the CF and TS envelopes to determine the CF and TS amplifica- 
uon factors at transition. These amplification factors are summarized in table A-l, and figure A23 
shows the traces of the amplification factors for a distance of 2.5% chord on either side of the 
experimentally determined transition location. These traces show the sensitivity of the correlation 
of experimental transition with analytically determined disturbance growth. 

As a supplement to the F-14 VSTFE cases, a few cases from earlier NLF flight tests (F-lll 
(ref. 8) and 757 (ref. 9)) were reanalyzed by the USS. The reanalyses of these data were consid- 
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ered important because these two sets of data had been used previously to define a transition 
criterion, but the boundary layer analysis and stability analysis procedures of the USS represent an 
improvement over that used in those previous studies. The cases chosen tend to define the enve- 
lope of the F-111/757NLF transition data band. Pressure distributions and disturbance growth 
curves for those cases are shown in figures A24 through A4 1 . Tables A-2 and A-3 summarize each 
case with their updated results. Figure A42 shows the movement of the transition point N-factors 
for the F-lll and 757 cases due to the change in analysis methods. 

Transition N-factors from the present study and those recalculated from previous data are pres- 
ented together in figure A43. The mean and estimate of standard deviation using y^^ 2 _j_ 2 

are also shown in figure A43. The scatter in the transition N-factors is large, even for the VSTFE 
points, for which the experimental data are most complete and accurate, indicating the need for 
further improvement in transition correlation. The USS should be a useful tool for further re- 
search. 


76 



Table A-1. — Summary of F-14 VSTFE Cleanup Glove Cases 
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200.8 



Hot-film locations 

(Only one row was installed at any one time) 


Figure A1. — Plan View of the Glove Showing Pressure Isobars, Hot Film Locations, and 
Static Pressure Geometry 
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(a) Pressure Distribution 


- CF (a\ s : 425 to 1,181 I/ft) 
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Case VI 

O Upper surface 
O Transition at x/c = 0 
O LE sweep = 35 deg 
O Re c = 10.79 million 
O Mach No. = 0.703 
O a = 1.1 deg 
O Flight 13, PT 1C 





(a) Pressure Distribution 



(b) Compressible Stability 


Figure A4. — Pressure Distribution and Compressible Stability Results for Case V2 




.0 


Input Points 


(a) Pressure Distribution 

CF {a * f# : 452 to 1,816 I/ft) . 

(Frequency: 0 Hz) 

TS (Frequency: 3,971 to 18,669 Hz) 

(Wave angle: 50 deg) 

/ ° 

Case V3 
Upper surface 

/ ° 

Transition at x/c = 0 


LE sweep = 15 deg 

y ° 

Re c = 24.77 million 

^y^ o 

Mach No. = 0.802 

y 00 ^ ° 

a = 0.2 deg 

o 

Flight 17, PT 31 B 


(b) Compressible Stability 

Figure A5. — Pressure Distribution and Compressible Stability Results for Case V3 
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Figure A6. 



Compressible Stability Results for Case V3 Including Maximum Amplification 
Method 
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(a) Pressure Distribution 
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Upper surface 
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LE sweep = 20 deg 
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Flight 18, PT 17B 


(b) Compressible Stability 

Figure A7. — Pressure Distribution and Compressible Stability Results for Case V4 





Figure A8. Pressure Distribution and Compressible Stability Results for Case V5 
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(a) Pressure Distribution 

CF (a\ s : 638 to 2,350 I/ft) 

(Frequency: 0 Hz) 
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(a) Pressure Distribution 

- CF (a*,,.: 432 to 1,823 I/ft) 
(Frequency: 0 Hz) 

■ TS (Frequency: 2,314 to 21,390 Hz) 
(Wave angle: 55 deg) 


Case V7 




O Upper surface 
O Transition at x/c = 0 
O LE sweep = 31 deg 
O Re c = 12.14 million 
O Mach No. = 0.802 
O a = 2.4 deg 
O Flight 13, PT 10A 
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(b) Compressible Stability 

Figure All. — Pressure Distribution and Compressible Stability Results for Case VQ 
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Figure A12. — Compressible Stability Results for Case V8 Including the Maximum 
Amplification Method 


.0 


Input Points 


(a) Pressure Distribution 



CF 

TS 

(a\ s : 498 to 1,654 I/ft) 
(Frequency: 0 Hz) 

(Frequency: 1,116 to 17,486 Hz) 


(Wave angle: 30 deg) 

/ 

O 

Case V9 
Upper surface 


/ 

O 

Transition at x/c = 0 


/ 

o 

LE sweep = 30 deg 


/ 

o 

Re c = 10.53 million 


/ 

o 

Mach No. = 0.704 


/ 

o 

a = 1 .4 deg 


./ 

o 

Flight 13, PT 2C 


(b) Compressible Stability 

Figure A13. — Pressure Distribution and Compressible Stability Results for Case V9 





(a) Pressure Distribution 



(b) Compressible Stability 


Figure A 14. — Pressure Distribution and Compressible Stability Results for Case W1 
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(a) Pressure Distribution 
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Figure A15. — Pressure Distribution and Compressible Stability Results for Case W2 
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(a) Pressure Distribution 
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Figure A1 6. — 


(b) Compressible Stability 

Pressure Distribution and Compressible Stability Results for Case W3 
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Case W4 

O Upper surface 
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(a) Pressure Distribution 

■ CF («' r# : 599 to 1,8541/ft) 

(Frequency: 0 Hz) 

■ TS (Frequency: 7,086 to 14,157 Hz) 
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Case W5 


Upper surface 
Transition at x/c = 0 
LE sweep = 25 deg 
Re c = 21.06 million 
Mach No. = 0.800 
a = 1.5 deg 
Flight 18, PT 18D 


(b) Compressible Stability 

Figure AT 9. — Pressure Distribution and Compressible Stability Results for Case W5 




Figure A 20. — Compressible Stability Results for Case W5 Including the Maximum 
Amplification Method 
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x/c 

(b) Compressible Stability 

Figure A21. — Pressure Distribution and Compressible Stability Results for Case W6 
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(a) Pressure Distribution 
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(b) Compressible Stability 

Pressure Distribution and Compressible Stability Results for the 757 NLF 
Glove Case 2 


BE 




N 


Figure A 25. 
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Compressible Stability Results for 757 NLF Glove Case 2 Including the 
Maximum Amplification Method 
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Figure A26. Pressure Distribution and Compressible Stability Results for the 757 NLF 

Glove Case 9 
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Figure A27. — Compressible Stability Results for 757 NLF Glove Case 9 Including the 
Maximum Amplification Method 
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Figure A28. — Pressure Distribution and Compressible Stability Results for the 757 NLF 
Glove Case 1 1 
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Figure A29. — Compressible Stability Results for 757 NLF Glove Case 1 1 Including the 
Maximum Amplification Method 
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(b) Compressible Stability 


Figure A30. — Pressure Distribution and Compressible Stability Results for the 757 NLF 
Glove Case 79 
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— Compressible Stability Results for 757 NLF Glove Case 19 Including 
the Maximum Amplification Method 
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Figure A32. — Pressure Distribution and Compressible Stability Results for F-111 
NLF Glove Case 12 
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Figure A33. 



Compressible Stability Results for F-111 NLF Glove Case 12 Including 
the Maximum Amplification Method 
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(b) Compressible Stability 

Figure A36 — Pressure Distribution and Compressible Stability Results for F-111 
NLF Glove Case 20 
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Figure A37. 


— Compressible Stability Results for F-111 NLF Glove Case 
the Maximum Amplification Method 
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Input Points 


(a) Pressure Distribution 




(«' r$ : 597 to 1,7101/ft) 
(Frequency: 0 Hz) 

(Frequency: 4,863 to 18,145 Hz) 
(Wave angle: 40 deg) 


F-1 1 1 NLF Case 22 
O Upper surface 
O Transition at x/c = 0.1 
O LE sweep = 25.2 deg 
O Re c = 29.69 million 
O Mach No. = 0.848 
O a = 4.28 deg 
O Flight 158, run 23 


(b) Compressible Stability 

Figure A38. — Pressure Distribution and Compressible Stability Results for F-1 11 
NLF Glove Case 22 






Figure A39. — Compressible Stability Results for F-111 NLF Glove Case 22 Including 
the Maximum Amplification Method 
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(b) Compressible Stability 

Figure A 40. — Pressure Distribution and Compressible Stability Results for F-111 
NLF Glove Case 24 
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Figure A41. — Compressible Stability Results for F-111 NLF Glove Case 24 Including 
the Maximum Amplification Method 


Table A-2. — Summary of 757 NLF Glove Cases 



05 


I s - 

O) 

in 

CD 

C D 

cvi 

CO 

CM 

CM 

C T> 


o 

r- 

CO 

oJ 

T— 

in 



CO 

CM 


CD 

CO 

05 

o> 


LO 

cvi 

LD 

CM 

I s - 

CO 

CO 





CM 

o 

1 

I s - 

CO 

CD 

T“ 

I s - 



CO 

r*- 

05 

05 

o> 

05 

CD 

CD 

CO 

CD 





CO 


o- 

o 

in 

o 

O 

00 


CD 

00 


CO 

CO 

05 

o 

CO 


O) 




CO 

m 

CM 

00 

CM 

CO 

T“" 

o 

CM 

r- 


rt 

LO 


in 

in 

CD 

m 

m 

in 

05 

o 

CO 

o 

I s - 


T~ 


CM 


CO 

CO 

m 

O 

o 

CO 

o 

cvi 

CO 

00 

CM 

CO 

CD 

CO 

00 

CO 

05 

CD 

in 



o 

I s - 

in 

T — 

CM 

00 


CD 

LO 

in 

CM 

05 

in 

CO 

Tf 

CO 


o 

o 


o 

V— 

in 

in 

o 

o 

o 

o 

o 






CM 

m 


05 


CM 

in 


T 



in 

CM 

in 

CO 

05 

co 

CO 

CO 

00 

CO 

csi 

CM 

c\i 

cvi 

CM 

05 

05 

I s - 

CO 

I s - 

CM 


CM 

O' 

rr 

00 

CO 

CO 

CO 

00 


CM CD 05 


O) 

O 


CM 

CM 

in 

05 

CD 

ih 

LO 



V“ 

CM 

CM 


< 

in 

CO 

CO 

r^ 

CO 


CM 



r- 

r- 

CO 

in 

CD 



in 

in 


CM 

CD 

o 

CM 




CM 

CM 

CM 















14 r 


12 h 



□ 757 NLF 
A F-111 NLF 


Nts 


6 

4 

2 

0 


20 


4P& 


12 



10 


Movement of transition N-factors from 
the previously documented values to 
those obtained by the USS code for 
the noted analysis case numbers 
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Figure A42. — Change in Transition N-Factors Using USS 



Figure A43. — Transition Mean Line and Standard Deviation From Mean 
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